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SUMMARY 

A FORTRAN program that can be used to analyze the elastic scattering of spin 
zero particles from an atomic nucleus is presented. The nucleus is represented 
by the nuclear optical model in the form of a Saxon-Woods vol\mie potential plus 
a surface Interaction potential whose radial dependence is in the form of the de- 
rivative of a Saxon-Woods form factor. There are available four independent pa- 
rameters of each of the two potentials; V and W, the strengths of the real and 
imaginary parts of the potential, R, the radius of the potential, and a, the 
diffuseness parameter. Provision is made for varying these parameters over a 
chosen grid in the parameter space or for searching for a best fit to experimen- 
tal cross section data by minimizing the chi- squared deviation as a function of 
chosen parameters. Since most of the computation time is expended in integrating 
the radial wave equation, the method of integration chosen is about twice as fast 
as the commonly used Runge-Kutta method. 


INTRODUCTION 

This report describes a FORTRAN program, ELSA, developed at the Lewis Re- 
search Center for the analysis of elastic scattering of spin zero particles 
against atomic nuclei, which are represented by the nuclear optical model. ELSA 
was developed to analyze experimental scattering cross section data in terms of 
the optical-model parameters, and, for this purpose, includes a search procedure 
for obtaining good fits to experimental data. The calculations are repeated many 
times in the search procedure, so the program was designed to use a method of 
integration (ref. l) that is about twice as fast as the commonly used Riznge- 
Kutta method (ref. 2). Since most of the computation time is expended in the in- 
tegration of the radial wave equation, this method leads to a relatively fast 
program. 

The program calculates the differential elastic cross section ^ 

function of the center of mass scattering angle G for spin zero particles with 
arbitrary mass, charge, and nonrelatlvistic energy scattered from a nucleus at 
rest with arbitrary mass and charge. The incident and target particles are taken 
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to have an interaction in the form of a tvo-hody potential, vhich is taken to he 
a complex nuclear potential containing a volume term of the Saxon-Woods (ref. 3) 
type, and a surface term whose radial dependence is in the form of the derivative 
of the Saxon-Woods foim factor. The spin of the target nucleus is neglected in 
accordance with the usual formulation of the nuclear optical model. If the in- 
cident particle is charged, the interaction includes a coulomh potential between 
an incident point charge and a target sphere of constant charge density. 

The potential parameters available are V and W, the real and imaginary 
constants that fix the depth of the potential well; R, the radius of the poten- 
tial well; and a, the constant that controls the diffuseness of the surface of 
the well. Each of these parameters can be specified independently for the volume 
and surface potentials. The parameters may be varied by fixed amounts by means 
of a grid variation to cover a specified set of points in the eight-dimensional 
parameter space. It is also possible to search the parameter space for a set of 

parameters that minimizes (chi-squared deviation) computed from the experi- 
mental differential scattering cross section, calculated 


MATHEMATICAL FORMULATION 
General Scattering Formulas 

First, the basic equations of scattering for charged, spin zero particles 
will be reviewed. The interaction between the particles is taken to be of the 
form 


V = Vj^(r) + V^(r) 


( 1 ) 


where is the complex nuclear potential and is the coulomb potential, 
both depending only on the distance r between the incident particle and the 
target particle. The Schrodinger equation for this case is 


- g + V(r)jf = 


Ef 


( 2 ) 


■where 


fflp + niQi 


(3) 


E 


m^ 

nip -t- 


(4) 


and lOm are the mass of the incident and target particles, respectively, and 
®lah laboratory energy of the incident particle. 

The wave function for the incident particle must be a solution at large val- 
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ues of r to the Schrodinger equation with VTr(r) =0 and -with 




4jt€or 


( 5 ) 


vhere is the charge of the incident particle and Zipe the charge of the 

target. For convenience, the parameter r) is introduced: 


h^k47t€Q 


( 6 ) 


where k is the wave number in the center of mass system. The Incident wave 
function normalized to one Incident particle per unit time per unit area is de- 
noted hy it is also a solution to 

(7) 

In this case, 

= -4= r(l + e^^^F[lri, 1, ik(r - z)] (8) 

where v is the velocity of the particle, the z coordinate axis is taken as 
the direction of the incident particle, and F is the confluent hyper geometric 
function. 


The asymptotic form of can he written as 




_ 1 J i[kz-T] In k(r-z)] 




^c( l(kr-Ti 2kr) 


+ e 

r 


^ - A 

The Rutherford scattering amplitude here is denoted hy f^3(0) and is 

f- -n ln[sin2(e/2)]_2lCT^ 

2k sin^ I 


>} 


(9) 


( 10 ) 


The zero angular momentum coulomb phase shift is obtained by setting 

the orbital angular momentum quantum number L equal to zero in the usual cou- 
lomb phase shift formula: 

0 L = arg r(L + 1 + iri) ( 11 ) 

The partial wave expansion of is 
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( 12 ) 


1 \ ^ T ioT 0) 

^ (2L + l)i^e %L(ii,kr) 

L=0 

where is the couloinb phase shifty Fj^(T)^kr) is the regular coulomb function, 

and Pj^(cos 0) are the Legendre polynomials. 

The total wave fimction ^tot written as the sum of the incident 

wave and a scattered wave, where the scattered wave includes only deviations from 
pure Rutherford scattering and interference terms. This can be written in the 
partial wave expansion: 


I \ ^ T io-T Pt(cos e) 

*tot - ’■= * *s=att 

L=0 

The radial wave function must satisfy the condition that only the out- 

going wave is modified by the nuclear interaction, which means that the asympto- 
tic expression is 


'^L = i’L('nAr) 



+ 



(14) 


•where is a complex constant and is the irregular coulomb function. 


When the asymptotic forms for the regular and irregular coulomb functions 
are introduced. 


Fj^ s sin^kr - t) In 2kr “ ^ + o'l) 
= cos^kr - T) In 2kr ~ ^ + <^ 1 ,^ 


(15) 

(16) 


the radial 'wave function becomes 


/ Lit \ ®L i [kr-T] In 2kr-(Lrt/2)+a "1 

'Fl = sin(kr - t] In 2kr - ~ + crhj ® ^ ^ 


It is possible to express another manner in the asymptotic region, name- 
ly, as a regular coulomb vave whose phase has been shifted by an amount gi-ven by 
the nuclear phase shift This expression is 


'Hr 


= sin/kr 


In 2kr 






(18) 
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where is a complex constant and 5 l the complex nuclear phase shift, A 

comparison of the coefficients of the incoming and outgoing waves in equations 
( 17 ) and (18) yields 


- e 


215t 


At = e 


±br 


(19) 


J 


Thus in the asymptotic region can he expressed as 


^ sinikr - 


In 2kr 


Lit ^ \ 

“ 2 


sin bjB 


i|kr-Ti 1^ 2kr-(L5t/2)+a]^+hJ 


(20) 


The asymptotic form of the wave function can he written in terms of the asympto- 
tic form of the incident wave and the coefficient Bj^ as follows: 


1 f i[kr-T] In k(r-z)] 


^ L ik(r - z)Jj 


^ gl(kr-ii In 2kr) 
-y/v n 


00 

fc(®) + ^(2L+ ^ P^Ccos 6) 

L=0 


ZIot. 


( 21 ) 


Ikr / 

The scattered wave is the portion containing the e ^*^/r term, and the dif- 
ferential cross section is given hy 




( 22 ) 


where 


%(e) = i^(2L + l)e 
L=0 

and f(,(9) is given hy equation (lO), 

The pure coulomb cross section is 


2i0T 

^ il ^l(cos e) 


(23) 


a^O) = \f^{e)\ 


(24) 
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Nuclear and Coulomb Potentials 
The nuclear potential is taken to he 


'r„(r) . - 


+ iW, 


1 + exp 


r - Ri 


(Vg + 


(l + exp 



( 25 ) 


The numbers Vg and Wg are positive for an attractive potential. 

The coulomb potential is taken to be that of a sphere of radius having 

a uniform charge density: 


„ 1 L r^ 1 

''c - 2 I " r2 , 


4rt£or 


r < Rr 


r > R, 


(26) 


PROGRAM DESCRIPTIOR 
General Description 

The program is written in FORTRAN II language and^ as listed in appendix 
is arranged to operate with the Lewis Research Center 7090 monitor system. The 
function of the monitor is to coordinate compiler and assembler processing and to 
provide a means for initiating execution of object programs. The monitor system 
also provide an error package whose function is to detect and control certain 
types of errors in FORTRAN programs. The error control cards used with ELSA are 
designated *MAJ0R UNDER (NO PRINT) and *ZER0 DIVIDE, These correct a floating- 
point major underflow and a division by zero by setting the result to zero. The 
library tape functions used are the standard FORTRAN II functions^ with the ex- 
ception of a subroutine called PLOTMT (X^Y^K^P). This subroutine is part of the 
library tape of the Lewis Monitor System (ref. 4) and is available for general 
use. Where this subroutine is not on the library tape, the program can be run by 
changing statement 715 (card 7-073) to read ”RRTURN, " deleting cards 7-074 and 
7-075 in subroutine SEVEN, and deleting subroutine PL0T7. The last input data 
card, with the data for the horizontal and vertical scales of the plot, will not 
be used if this change is made. 

ELSA consists of a main routine and 10 subroutines. The main routine calls 
the subroutines in the proper sequence and also provides for a sequential varia- 
tion of the parameters Vj^, Rj^, and a^_ (i =1,2). The subroutines and the 

calculations each perfoimis are listed as follows: 


6 


Subroutine ONE 
Subroutine TWO 
Subroutine THREE 
Subroutine FOUR 

Subroutine FIVE 
Subroutine SIX 

Subroutine SEVEN 
Subroutine PLOT? 
Subroutine PLOTMY* 

Subroutine EIGHT 


Reads in data cards, performs a reduction of input data 
Computes coulomb phase shifts and coulomb wave fionctions 
Tabulates coulomb potential and r^ 

Tabulates optical potential form factors, obtains power- 
series coefficients 

Tabulates effective kinetic energy 

Computes radial wave function by integrating Schrodinger 
equation, then computes scattering amplitudes 

Computes cross section X^, calls PLOT? 

Sets up PLOTMY 

Furnishes a plot of a against 6 , and a, , against 

ex ex un 

®th 

Provides a search on any combination of the parameters Vj^, 
Wj_, Rj^, and aj^ (i =1,2) by minimizing X^ 


Details of Specific Subroutines 

Subroutine ONE . - This routine reads in the data cards and performs a reduc- 
tion of the input data to the forms used by the program. It determines the size 
of the increment (= 6, FORTRAN name ADD) used in the integration. This is 

done by setting 5 = Sq/f and rounding the result to three figures. The dimen- 
sionless input parameter 6^ = DEL is approximately the size of the integration 
increment in terms of the dimensionless variable kr. The subroutine also com- 
putes normalization factors for the radial wave fimctions. The input data used 
are listed in the section INPUT DATA. 

Subroutine TWO . - The coulomb phase shifts are computed by using Stirling's 
formula to obtain a value for a^Q from equation (ll) and then by \ising the re- 
currence relation 

^L-1 “ ^L ” 'tan"^ ^ = AS(L) (2?) 

to obtain the values for the lower values of L. The designation listed to the 
right of the (=) sign is the FORTRAN name of the variable. 

The regular and irregular coulomb wave functions are computed at two 

values of r, namely, r^^^^ and r^^^ - 26 (ARM and ARM -2 ADD). In order to 

furnish starting values for the use of the recurrence relation, the values of the 
coulomb wave functions for L = 0 and L = 1 (Fq, Gq_) are computed at these 


? 



points by using the asymptotic relations (ref. 5) 


cos sin 0^ = FF(N, L + 1) 

Gl = cos 0j^ - sin 0^ = GG(N,L + 1) 
vhere N is equal to 1 or 2 depending on the radial position, and 

ttT I 

0j^ = kr - Ti In 2kr ~ ~ 

The variables and used to calculate the coulomb wave functions are 

given by the series 


( 28 ) 


(29) 


Ul + iVj, = (1 + C]_^L *^1,L^2,L '^1, L^2, L^3, L 


+ . . 


■'J,L - 


_ (in - L + J - i)(in + L + J) 


i2krJ 


(30) 

(31) 


The maximum number of terms in the asymptotic series is set by the input control 
Ml, The series is terminated vhen the magnitude of any term exceeds that of the 
previous term, or when the ratio of the magnitude of any term to the magnitude of 
the sum of the series becomes less than the input number BBS. If neither of 
these terminations occurs and the number of terms actually reaches Ml, the value 
of AEM is increased to AEM + 6 X M2 x ADD, and a new series is computed. The 
values of Fq, Fq_, Gq, and are checked by computing the Wronskian, and if 

they are not sufficiently accurate, the value of AEM is increased to 
AEM + 6 X M2 X ADD, and the coulomb functions are computed at this new point. 

The values of Fj^ and Gj^ are then obtained from the recurrence relation 

lV(l+ 1)2 + ULfl = (2L + 1 )[ti + - (L + 1) y/iFT^ Ul_i 

(32) 

where 

Ul = Fl or Gl 

The recurrence relations can be used to obtain accurate values of G^j^ by re- 
curring upward from Gq and G^^. Since the Fj^ cannot be computed accurately 
by recurring upward, arbitrary starting values are assigned to I^lEMMO 
Flemmp ■where LEM is the maximum value of L to be used in the program and is 
determined from input parameters. The recurrence relations are used to generate 
the F;|^ downward,, which are then normalized to the value of Fq computed by the 
asymptotic series. Each F^ is then checked to see if the Wronskian satisfies 
the following criterion: 
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(33) 





< 10 "^ 


If any of the fail this test, the starting value of L for the dovnvard re- 
currence on is increased hy 10, and the process is repeated. If the start- 

ing value of L for the dovn^ard recurrence becomes greater than L = 100, the 
program value of APM is increased and the following statement written out: KE- 

CUEREINCE RELATIONS FOR COULOMB WAVE FUNCTIONS HAVE FAILED, INCREASE RMAX, If the 
value of ARM is ever increased to the point where the number of storage locations 
allocated to and r^ in subroutine THREE will be exceeded, the program pro- 
ceeds to the next set of data, and the following statement is written out: ARM 

EXCEEDS AVAILABLE STORAGE FOR POTENTIALS, PROCEED TO NEiXT CASE. 

Subroutine THREE . - The values of the coulomb potential V^ and of r^ are 
tabulated at all mesh points. The coulomb potential is added to the constant 
kinetic energy term and stored in the forms 




r < R, 


12 V 


1,2 . . VC(I) 


r > R, 


(34) 


Subroutine FOUR . - The Saxon-Woods form factor 

1 


fl(r) = 


1 + exp 


^ 


= SS(I) 


(35a) 


is tabulated at each mesh point. The surface interaction form factor. 


4 exp 


fp(r) = 


- ^2 

_ ^2 . 


{• * 


= TT(I) 


(35b) 


is also tabulated at each mesh point where it is significant. 

To start the integration of the radial wave function in subroutine SIX, a 
power-series solution is used to get initial values of the wave function. This 
solution uses a power-series expansion for the nuclear potential at a starting 
point that must be inside the nuclear potential well, since the series converges 
only for r < Rq_. Since the starting value of r is usually small, the surface 
potential is neglected in obtaining the initial value of the wave function. The 
coefficients for the power-series expansion of the Saxon-Woods form factor are 
computed by using 
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1 


(36) 


1 + exp 


"r - Ri“ 

^-1 


00 

= y r 


m=0 


The coefficients can he expressed in terms of the form factor evaluated 

at r = 0. The exponential evaluated at r = 0 is denoted hjr 


exp 


[i]= 


E 


(37) 


The are then given hy 


rm = 


=5^ (t^) 


J=1 


M 


The coefficients Aj satisfy the recurrence relations 

Aj = JAj‘^ + Aj”^ = AC(J) 

with 


.M 


^1 = 


1 = A 


M 


(38) 


(39) 


and 


aM =0 


J > M 


= 0 


J < 1 


Subroutine FIVE . - The nuclear potential contribution to the effective ki- 
netic energy term in the Schrbdinger equation is obtained and added to the cou- 
lomb potential. This is stored in the following forms: 


VC(I) + 


ft2 12 


2ja ^ 

fi2 12 


[Vifi(r) + V2f2(r)] 
[wifi(r) + W2f2(r)] 


= W(I) 


= WW(I) 


(40a) 

(40b) 


where VC(l) is computed in subroutine THREE, and fi(r) and f2(r) are computed 
in subroutine POUR. The power-series coefficients for the nuclear form factor 
are combined with the strength of the potential and stored as 
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(41a) 


DD(1,M + 1) 


DD(2,M+ 1) (411d) 

Subroutine SIX , - The radial vave functions are computed for each value of 
L from zero to integrating the Schrodlnger equation, starting near the 

origin and proceeding out to a point *vrell heyond the range of the nuclear poten- 
tial. The starting values for the “wave function are obtained from a poTver- series 
solution to the radial wave equation at the starting value of r. The differ- 
ential equation for the radial vave fiinction is 


ft ''1% = 

2m- tt - 

^ - 


_dr^ 


, 3.2 _ k2 XL = 0 


E 




(42) 


where 


Xl = = C 


00 


(43) 


n=0 


in vhich is the normalization constant given in equation (48). The effec- 

tive kinetic energy terms are expanded in powers of r as follows: 

^ m=0 

The following recurrence relations are found for the 0^: 


= n(n+ 2 L V l)|-^^®n-2 + ^ 4-2 “ § 4-4 ~ 

L c 



('m'^n-2-m “ ''^m®n-2-ii^ 


( 45a) 


{- 


= n(n + 2L“rTy|"^"4-2 + ®n-2 ~ ^ 


n-2 -j 

n-4 " ^|^^('m®n-2-m ^ ^m‘*’n-2-m) 
m=0 


The quantities in equations (45) appear in the program as 

$n = 4 i$n = BB(1,K + l) + 1BB(2,N + l) 


(45b) 


(46) 
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and 


, v2 4 v2 

Vm + iv" = + ^ ViTm + ^ = DD(l,M + 1) + iDD(2,M + 1) (47) 

The values for are tahen to he 

•I-Q = 1.0 ^0 = ® 


The starting values for the vave function are normalized to the coulomh 
function normalization so that 





a/i^ + 

1(2 . 1+1) • * * L(2 • L + 1) 


(kr)^^$ 


Ct (kr)^^$ 

(48) 


The integration method is that of successive extrapolation; that is, the 
value of Xl at any given point is extrapolated from the values at two previous 

points and the difference in r. The relations are obtained from Taylor series 
expansions of XL(r±5) about the point r and from use of the differential equa- 
tion to evaluate the second derivative appearing in the expansion (ref. l). The 
complex function Xj, is denoted as 


Xt = Xt + i 




(49) 


The equation to be solved can be written with the subscript L temporarily 
suppressed in the form 


— ^ (r) + iT)"(r) (X"*" + iX“) = 0 


(50) 


If 


q(r) = 1 + 


12 


(r) 


P(r) 



n”(^) 


the solutions are 

y±f-r + 5) = X~(r)q(r + 5) ± X~^(r)p(r + S) 
q^(r + 5) + p^(r + 6) 


(51a) 

(51b) 


( 52 ) 


where 
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By comparing equation (50) to the previous form of the radial equation (eq. (42)), 
it is seen that 


V(r) = + ^|Re Vjj(r) -V(,(r)J - 


L(L + 1) 


(54a) 


T)“(r) = + -^ Im %(r) 


The integration is carried out to the place vhere r = r^^^, 
of the wave function at the last two mesh points are stored. 


(54h) 

and the values 


The value of the scattering amplitude is obtained by evaluating equa- 

tion (14) at each of the last two mesh points. _^en equation (14) is written ex- 
plicitly with an arbitrary normalizing factor the equation becomes 


Cl(^L + ( 21 - j (^L ^ 


(55) 


Using this relation at two values of r enables elimination of and solution 

for Bj^. For a specified value of L, this solution is 


+ IB" 


2i[ 

K * 

>^2 - 

+ iXg 

;)^i] 


(4 + 

)(Gi + iFj 

-> - ( 

Xi + iX£^ 

|( G2 

+ iPg) 


= BP + IBM 


(56) 


where the subscripts refer to the last mesh point (2) and the previous mesh 
point ( 1 ) . 


Subroutine SEVEN . - This subroutine computes the differential elastic cross 
sections as given by equations (22) and (23) for each of the experimental angles. 
A value of x2 is computed as follows: 


^ex,N 

^ ^ ®ex, i ) 

®=®ex, i 


<^th(®ex, i) " ‘^ex(®( 

r 


ejc^ex, i 


■^ex. 




(57) 


where W(0 ^ ., ) is an assigned weight factor^ is the calculated cross sec- 

0X^ 1 II 

tion, and Uex is the experimental cross section at each experimental angle 
®ex* 

If a search over the optical model parameters is desired, subroutine EIGHT 
is entered. At the end of the search, or if no search is used, the differential 


elastic cross section is computed at vhere 6^-^ is a calciilated set of 

evenly spaced angles if KDCAL ^ 0, and otherwise is equal to if 

KDCAL = 0, The pure coulomb differential cross section, as given hy equation 
(24), is computed at each angle and the angles, elastic cross sections, 

coulomh cross sections, and ratio of elastic to coulomb cross sections are 
printed out. 

Subroutine EIGHT . - The search routine provides a variation of any of the 
parameters Vj^, Wj^, aj_, and (i =1,2), which are denoted by aj in the fol- 
lowing discussion. First is changed to a]_ + daj_, where da^ is input 

data. The subroutine returns to the optical model program, and a new set of 

p p 

02 ( 0 ) and a value of Xg are calculated. The value of Xg is compared 'with 
X^^ and if it is smaller, the variable is changed to + 2daQ_ and the 

process repeated, (If > X^, the variation is made in the opposite direction, 
that is, Oq is changed to a-i - dx-i . ) The variations are continued in the same 
direction until the value of the current X is larger than the previous X . 

p 

The value of a]_ for which X is a minimum is then interpolated by fit- 
ting a parabola to the last three values of X^(a]_). 

With the value of a]_ fixed at the interpolated value of (and keeping 
the six other aj fixed at the initial values), the value of cx-g changed to 
Og + Values of and X are computed, and the same process as de- 

scribed for cc^ is carried out. The procedure is repeated for each of the de- 
sired ttj. At the end of one of these circuits, the process can be repeated if 
desired with the same da^ or with each daj reduced by a constant factor. 

The search can be repeated as many times as desired and can be terminated by 
limiting the number of circuits, or when X^ fails to be reduced by some prede- 
termined input factor. When the search is terminated, the program returns to 
subroutine SEVEN, which prints out (and plots) the final cross section. 

Subroutine FLOTMY (X,Y,K,P ) .(ref^.^* - This routine is a general plotting 
subroutine that is part of the Library Tape of the Lewis Monitor System. Its use 
here is to furnish plots of log^o crex against 0^^ and log^Q <^th against 
©^j^. In the PLOTMI subroutine, the variable X is plotted down the page, and Y 

is plotted across the pagej K is an array controlling the possible options of 
PLOTMY, and P is an array controlling the vertical and horizontal scales. 

The subroutine PLOTMY is called by subroutine PLOT?, which prepares the 
cross-section data in a foimi suitable for PLOTMY. 

Subroutine PLOT?. - This subroutine arranges to put the values of 
-log-]_Q CTgx(®) ^^<3. ^og^o ‘^th*^®) array X(l) and the values of 9^^ and 

©.j-j^ into the array Y(l). The proper values of K(l) are set to plot two curves 
with the X and Y scales specified by input controls P(l). The input con- 
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trols P(l) are read in "by subroutine PLOT?. If the subroutine PLOTMI is not 
available to the program, subroutine PLOT? should be removed along vith cards 
7-074 and 7-075. Card 7-073 should be changed to read 715 RETURN. The input 
data card “with the plot controls, P(l), will not be used if this change is made. 

The input controls must be integers. A convenient method of selecting the 
control numbers is first to choose the desired increment in the horizontal 

scale that will be equal to one printing position (there are 100 printing posi- 
tions horizontally on the page). P(ll) must be chosen to be an integer that is 
^ times a power of 10. P(9) is then the value such that 6 - P(9) gives the 

correct exponent of 10 in the equation 


P(ll) = 

The value of P(l0) is then obtained from the desired starting value (N) of Q 
by the equation 


P(10) = 

These values set up the horizontal scale. The values for P(8), P(6), and P(7) 
controlling the vertical scale are set up in a similar manner. There are 60 
line spaces per page on the vertical scale. But in contrast to the horizontal 
scale, which is limited to 100 points, the vertical scale has no such limitation 
and can be made to cover several pages continuously. 


INPUT DATA 
Description 

Units for the input data are as follows: 

Energy Mev 

Potential strength Mev 

Mass amu 

Length fermis 

Charge dimensionless multiples of electron charge 

Angle degrees 

Cross sections millibarns 
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I 

The data are input from tape 7, -which is prepared from IBM cards as follows: 

Card 1 FORMAT (13A6) title card - 78 Hollerith characters to identify the run 

Card 2 FORMAT (1515) control numbers (the subroutines in which these are used 

are in parenthesis) 

lOTl maximimi number of terms in asymptotic series for coulomb wave func- 

tions (2) 

KN2 number of triplets of mesh points to increase r^^^y if asymptotic 

series does not converge (2) 

KH3 0^ no surface interaction; yes s-urface interaction (4) 

]M4 maximum number of terms in form-factor power series (4) 

1TN5 0, no search; 1, yes search using subroutine EIGHT (7^8) 

KN6 number of search circuits before search interval is reduced (8) 

NN7 maximum number of circuits in search (8) 

KN8 maximum number of terms in wave-function power series (6) 

NTl number of angles for final calculated ang-ular distribution (7) 

HTT n-umber of experimental angles (7) 

INKODE 0, do not read in experimental data; 1, read in experimental data (l) 

KDPPLT 0^ no plot is given; yes plot 
KDGRB 0, no rough grid; yes rough grid 

KDCAL 0, 0 ^-^ = ^ 6 ^^ 

Card 3 FORMAT (1515) control number (this card can be blank if no grid varia- 
tion is used) 

F71X mmiber of grid points for (MAIH) 

HWIX number of grid points for (MAIN) 

NRIX mmber of grid points for R;j_ (MAIN) 

NAIX number of grid points for a^ (MAIN) 

NV2X number of grid points for Vg (MAIN) 

NW2X number of grid points for Wg (MAIN) 
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mzx 


number of grid points for Rg (MA.IN) 

NA2X number of grid points for (MAIN) 

Card 4 FORMAT (8F10.0) iDasic data 

AVI \^l\’ ^^0*1 strength of Saxon ^^ell (5) 

AWl |^l|> strength of Saxon well (5) 

ARl R^, radius of Saxon well (4) 

AAl a^_, diffuseness of Saxon well (4) 

DVl grid or search increment of V^ (8 or MAIN) 

DWl grid or search increment of (8 or MAIN) 

DRl grid or search increment of R^^ (8 or MAIN) 

DAI grid or search increment of a^ (8 or MAIN) 

Card 5 FORMAT (8F10.0) basic data 

AV2 |^2l> strength of surface well (5) 

A ¥2 I ^2 1 ^ strength of surface well ( 5 ) 

AR2 Rg, radius of surface well ( 4) 

AA2 ag, diffuseness of surface well (4) 

DV2 grid or search increment of Y 2 (8 or MAIN) 

DW2 grid or search Increment of Wg (8 or MAIN) 

DR 2 grid or search increment of R^ (8 or MAIN) 

DA2 grid or search increment of a 2 (8 or MAIN) 

Card 6 FORMAT ( 8 FIO.O) basic data 

AEE laboratory energy of incident particle (l) 

ARC radius of charge distribution (3) 

DEL desired mesh size 6 q as a fraction of kr (l) 

AER radius at which mesh size is doubled ( 6 ) 

AMP mass of incident particle (l) 
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charge of incident particle (l) 

AMT mass of target particle (l) 

AZT charge of target particle (l) 

Card 7 FORMAT (SFIO.O) "basic data 

ATI initial angle for calculated angular distribution (?) 

DTI increment in angle for final angular distribution (?) 

Card 8 FORMAT (8E10.0) convergence parameters 


BBl) 

determine 

^max ^ + ^ 1®1 + 

^2 

( 1 ) 

BB2) 

determine 

^max ^ ^2 ^2®1 

®2 

(1) 

BBS 

determine 

^max > ®3 + 

( 1 ) 


BB4 

determine 

^ + ®4 

( 1 ) 



BBS convergence parameter for asymptotic series for coulomb wave func- 

tions (2) 

BBS form-factor cutoff (4) 

BB7 convergence parameter for wave-function power series (S) 

BB8 convergence parameter for scattering amplitudes (S) 

Card 9 FORMAT (8EIO.O) convergence parameters 
BB9 reduction factor for search interval (8) 

BIO convergence factor on decrease in by search program (8) 

Bll determine starting radius for integration from r = B -[-| L(5^/k); note 

Bll should always be greater than 0.4 (6) 

If INKODE = 0, there are no cards with experimental data to be read in. If 
lEKODE f 0, cards containing experimental data will be required. The nuimber of 
these cards is as follows: 

NTT even, there are NTT/2 cards 

NTT odd, there are (NTT + l)/2 cards 

Each card contains two angles, two cross sections, and two weight factors. 
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Cards 10, etc. FORMIAT (6F10.0) 

HH(1) first experimental angle (?) 

XX(l) first experimental cross section (?) 

WX(l) first veight factor (?) 

HH(2) second experimental angle (?) 

XX(2) second experimental cross section (?) 

■WX(2) second T^eight factor (?) 

If KDPLOT ^ 0, there is one more card: 

Card PLOT FOPMAT (SFIO.O) controls for PLOT? subroutine 

P(g) determines vertical scale factor, 10^“^^^^ 

P(?) determines vertical scale starting value; let M he hhe first 

integer greater than log^^Q a^gx^ then P(?) is the negative 
number such that |P(7)| =Mxlo6-P(S) 

P(8) determines desired increment AClog^o vertical scale 

equal to one line space, P(8) = A(log^Q a)l0^"^^^) 

P(9) determines horizontal scale factor, 10^“^^^) 

P(10) determines horizontal scale starting value; let N he desired 

starting value for 0, then P(10) = Nxl0^“^(^) 

P(ll) determines desired increment AQ in horizontal scale equal to 

one printing position, P(ll) = A8xl0^“^(^) 


Sample 

A sample data set of input data is listed in appendix B. There are nine 
cards corresponding to cards 1 to 9 as described in the previous section, 18 
cards containing the experimental angles, cross sections, and weight factors, and 
a final card with the PLOT controls. The values of the various controls numbers 
and convergence parameters, which appear on cards 2, 8, and 9, are t3rpical values 
that have been foimd satisfactory and that can be used as guides in choosing 
these parameters. 


SAMPLE OUTPUT DATA 

A sample set of output data from the program is listed following the sample 
input data listing (appendix B). Not all the possible output described appears 
in the sample case. 
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The title card is printed at the beginning of each case. Following this on 
the first page is a labeled print-out of the input data contained on cards 1 to 9. 
(Card 3 is printed here only if a parameter grid is being used. ) The experimen- 
tal angles, cross sections, and weight factors are printed on the second page. 

At the top of the third page some of the output from subroutine OKE appears. 
Here are printed values of the wave number k = AKKj the mesh size, Ar = ADD; the 
coulomb parameter, t) = AZZ; the normalization factor for the L = 0 radial wave 
function, ACC; the radius at which the mesh size is doubled, ARE; the matching 
or maximum radius, ARM; the value of the index of the doubling radius, IKR; the 
value of the index of the maximum radius, IKM; and the maximum L value to be 
used, LRM. 

The next output comes from subroutine TWO. If the recurrence relations give 
values of Pl Gl whose Wronskian is greater than 10“5, there is an output 

statement noting this and noting that ARM has been Increased. If the value of 
ARM is Increased beyond the storage assigned to the potentials, there will be an 
output statement noting this fact, and the program will proceed to the next set 
of data (if any). Since it is possible for ARM to be increased with no output 
statement, there is always a write- out of the final values of ARM and IRM at the 
completion of subroutine TWO. 

The final values of the coulomb radius ARC and the index IRC are written out 
from subroutine THREE after ARC has been adjusted to be at one of the mesh 
points. 


Next there is output from subroutine SIX. Each time the power series used 
to obtain a starting value for the radial wave function fails to converge, the 
starting value for r is decreased and a new power series is computed. The to- 
tal number of times this happens (for all L values) is printed out as MMM. If 
the value of r is ever decreased to a value too small for the mesh size used, 
the computation stops, the statement POWER SERIES FAILED TO CONVERGE IN SIX is 
printed out, and the current values of MMM and L are printed out. There is a 
provision for cutting off the number of partial waves used by the program when- 
ever the scattering amplitude becomes sufficiently small. Therefore, the maximum 
value of L that can be used is printed out as LRM, and the number of partial 
waves actually used is printed out as LMAX. 


If a search or a grid is being used, the values of the parameters and the 
value of y? = ERROR are printed at each point of the parameter space. When the 


search program is being used, the interpolated values of the yr 
eter being varied are printed whenever an interpolation is made. 


and the param- 


When a calculation with a single parameter set or with a search procedure is 
completed, the calculated values of the elastic cross section, the coulomb cross 
section, and the ratio of elastic to coulomb cross section, along with the cor- 
responding center of mass angles, are printed out. A plot of the cross sections 
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against angle is printed if desired. No print-out or plot of the cross section 
appears if a parameter grid is used. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, October 9, 1963 
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APPENDIX A 


PROGRAM LISTING 

C MAIN PROGRAM FOR *ELSA* M-000 

COMMON AAl »AA2*ACC*ADD.AKK,AKS,ARCfARM»ARR*ARl » AR2 » ASS t ATI ♦ AVI » C-01 

1 AV2*AW1,AW2,A2D»AZ7*A2D,A3D,A6D» C-02 

2 BE1»BB2,BB3,BB4,BB5*BB6,BB7»BB8»BB9»B10»B11»B12.BM,BP» C-03 

3 CC»CN»CS» C-04 

4 DAI ,DA2*DD»DDS»DR1 *DR2»DT1 »DV1 *DV2*DW1 ,DW2*D12. C-05 

5 ERR,FF»GG>HH»HS* C-06 

6 I I1*II2»I I3,II4*ICC»INK0DE,IRM,IRR»I3I , C-07 

7 KDCAL»KDGRD,KDPLT»KDDE2»KaDE6,KDDE7»KQDE8*K800» C -08 

8 LL*LLM*LRM» C-09 

9 NA1X*NA2X*NN1,NN2,NN3,NN4*NN5,NN6,NN7,NN8*NNN,NR1X»NR2X C-10 

COMMON NTN,NTT*NT1*NV1X*NV2X»NW1X*NW2X*N3N»PL* C-11 

1 QQ»RR*SN*SS,TT»VC»VV>WW*WX»XCf XN»XX C-12 

DIMENSION BM( 51) »BP(51 ) » D-01 

1 CC(2,100) »CN(90) »CS(51) » D-02 

2 DD(2,100) ,FF(2,51 ) »GG{2,51 ) ,HH(90) »HS(90) *PL(51 ) * D-03 

3 QQ( 50) »RR( 500) ,SS( 500) » 5N ( 51 ) * TT ( 500 ) fVC(500) *VV( 500) » D-04 

4 WW( 500) »WX(90) » D-05 

5 XC(90) »XN(90) »XX(90) D-06 

1 CALL ONE M-OOl 

2 CALL TWO M-002 

K0DE2=K0DE2 M-003 

GO TO (1,3)*K0DE2 M-004 

3 CALL THREE M-005 

IF (KDGRD) 10»4fl M-006 

10 TV1=AV1 M-007 

TW1=AW1 M-008 

TR1=AR1 M-009 

TA1=AA1 M-OlO 

TV2=AV2 M-Oll 

TW2=AW2 M-012 

TR2=AR2 M-013 

TA2=AA2 M-014 

12 DO 20 NA2=1»NA2X M-015 

IF (NA2-1) 31»30*3x M-016 

30 AA2=TA2 M-017 

GO TO 32 M-018 

31 AA2=AA2+DA2 M-019 

32 DO 20 NR2=1»NR2X M-020 

IF(NR2-1) 34»33,34 M-021 

33 AR2=TR2 M-022 

GO TO 35 M-023 

34 AR2=AR2+DR2 M-024 

35 DO 20 NW2=1»NW2X M-025 

IF (NW2-1) 37»36»37 M-026 

36 AW2=TW2 M-027 

GO TO 38 M-028 

37 AW2 = AW2+DVv'2 M-029 

38 DO 20 NV2=1»NV2X M-030 

IF (NV2-1) 40*39*40 M-031 

39 AV2=TV2 M-032 

GO TO 41 M-033 

40 AV2=AV2+DV2 M-034 

41 DO 20 NA1=1*NA1X M-035 

IF(NAl-l) 43*42*43 M-036 

42 AA1=TA1 M-037 
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GO TO 44 

43 AA1=AA1+DA1 

44 DO 20 NR1=1*NR1X 
IF (NRl-1) 46»45»46 

45 AR1=TR1 
GO TO 47 

46 AR1=AR1+DR1 

47 DO 20 NW1=1»NW1X 
IF (NWl-1) 49»48»49 

48 AW1=TW1 
GO TO 50 

49 AW1=AW1+DW1 

50 DO 20 NV1=1,NV1X 

IF tNVl-1) 52*51»52 

51 AV1=TV1 
GO TO 53 

52 AV1=AV1+DV1 

53 CONTINUE 

4 CALL FOUR 

5 CALL FIVE 

6 CALL SIX 
KDDE6=K0DE6 

GO TO (1,7) ,K0DE6 

7 CALL SEVEN 
KDDE7=K0DE7 

GO TO (1,8,20) ,K0DE7 
20 CONTINUE 
GO TO 1 

8 CALL EIGHT 
K0DE8=K0DE8 

GO TO (4,5), K0DE8 
END 


M-038 

M-039 

M-040 

M-041 

M*-042 

M-043 

M-044 

M-045 

M-046 

M-047 

M-048 

M-049 

M-050 

M-051 

M-052 

M-053 

M-054 

M-055 

M-056 

M-057 

M-058 

M-059 

M-060 

M-061 

M-062 

M-063 

M-064 

M-065 

M-066 

M-067 

M-068 

M-069 



SUBROUTINE ONE 1-000 

PROGRAM ONE DATA REDUCTION 1-001 

COMMON AA1»AA2,ACC,ADD»AKK*AKS»ARC»ARM,ARR*AR1*AR2»ASS,AT1 ,AV1 t C-Ol 

1 AV2 ,AW1 ,AW2*AZD»AZZ,A2D,A3D,A6D» C-02 

2 BBl ,BB2,BB3,BBA»BB5,BB6»BB7,BB8 ,BB9 ,B10 »B1 1 >B12 ,BM ,BP t C-03 

3 CC,CN,CS» C-04 

4 DAI ,DA2,DD»DDS»DR1 »DR2 >DT1 »DV1 ,DV2»DW1 ,DW2,D12» C“05 

5 ERR,FF,GG»HH*HS* C-06 

6 II1,II2»II3,II4*ICC»INKDDE»IRM,IRR*I3I , C-07 

7 KDCAL»KDGRD,KDPLTfK0DE2,K0DE6»K0DE7,)<0DE8,K8 00> C-08 

8 LL,LLM,LRM» C-09 

9 NA1X,NA2X»NN1,NN2,NN3*NN4,NN5,NN6,NM7,NN8,NNN,NR1X»NR2X C-10 

COMMON NTN,NTT,NT1,NV1X,NV2X,NW1X»NW2X »N3N*PL, C-11 

1 QQ,RR,SN*SS,TT,VC,VV»WW9WX»XC,XN*XX C-12 

dimension BM{51)*BP(51), D-01 

1 CC(2*100) ,CN(90) »CS{51) * D“02 

2 DD( 2,100) ,FF(2,51 ) *GG(2»51 ) ,HH(90) *HS(90) »PL (51 ) * D-03 

3 QQ( 50) »RR I 500 ) ,SS( 500) , SN ( 51 ) > TT ( 500 ) , VC ( 500 ) , VV ( 500 ) , D-04 

4 WW(500) ,WX(90) , D-05 

5 XC(90) ,XN(90) ,XX(90) D-06 

dimension TITLE( 13) D-07 

READ INPUT TAPE 7,101, ( T I TLE C I) , I =1 , 1 3 ) 1-002 

WRITE OUTPUT TAPE 6,102, ( T I TLE ( I ) , I = 1 , 1 3 ) 1-003 

READ INPUT TAPE 7,120, NN 1 ,NN2 ,NN3 ,NN4 ,NN5 ,NN6 , NN7 , NN8 , 1-004 

1 NT1,NTT,INK0DE,KDPLT,KDGRD,KDCAL 1-005 

READ INPUT TAPE 7,120, NVlX ,NW1X ,NR1X ,NA1X ,NV2X ,NW2X ,NR2X ,NA2X 1-006 

READ INPUT TAPE 7,116, AVI , AWl , ARl ,AA1 ,DV1 ,DW1 ,DR 1 , DAI , 1-007 

1 AV2,AW2,AR2,AA2,DV2,DW2,DR2,DA2, 1-008 

2 AEE,ARC,DEL,ARR,AMP,AZP,AMT,AZT, 1-009 

3 ATI, DTI 1-010 

READ INPUT TAPE 7,118, BBl , BB2 ,BB3 ,BB4 ,BB5 ,BB6 ,BB7 , BBS , 1-011 

1 BB9,B10,B11 1-012 

IF( INKDDE ) 150,151 » 150 1-013 

150 read input tape 7,125, ( HH ( I) ,XX ( I ) ,WX (I ) , 1 = 1 ,NTT ) 1-014 

151 WRITE OUTPUT TAPE 6,103 1-015 

WRITE OUTPUT TAPE 6 , 104 , AVI ,AWl ,AR1 ,AA1 1-016 

WRITE OUTPUT TAPE 6 , 105 ,DV1 ,DW1 ,DR1 ,DA1 1-017 

IF( KDGRD) 140,141 ,140 1-018 

140 WRITE OUTPUT TAPE 6 ,144,NV1X ,NW1X ,NR1X ,NA1X 1-019 

141 WRITE OUTPUT TAPE 6 , 106 , AV2 , AW2 ,AR2 , AA2 1-020 

WRITE OUTPUT TAPE 6 , 107 ,DV2 ,DW2 ,DR2 ,DA2 1-021 

IF( KDGRD) 142,143 >142 1-022 

142 WRITE OUTPUT TAPE 6 , 145 ,NV2X ,NW2X ,NR2X ,NA2X 1-023 

143 WRITE OUTPUT TAPE 6 , 108 ,AEE , ARC ,DEL , ARR 1-024 

WRITE OUTPUT TAPE 6 , 109 ,AMP , AZP ,AMT , AZT 1-025 

WRITE OUTPUT TAPE 6 , 1 10 , ATI ,DT1 ,NT1 ,NTT 1-026 

WRITE OUTPUT TAPE 6 , 1 1 1 ,BB1 ,BB2 , BB3 ,BB4 1-027 

WRITE OUTPUT TAPE 6 , 1 1 2 , BB5 , BB6 , BBT ,BB8 1-028 

WRITE OUTPUT TAPE 6 , 1 1 3 ,BB9 , B 10 , B1 1 1-029 

WRITE OUTPUT TAPE 6 , 114 ,NN1 ,NN2 ,NN3 ,NN4 1-030 

WRITE OUTPUT TAPE 6 , 11 5 ,NN5 ,NN6 ,NN7 ,NN8 1-031 

WRITE OUTPUT TAPE 6 , 160 , I NKODE ,KDPLT ,KDGRD ,KDCAL 1-032 

IF( I NKODE) 153,152,153 1-033 

152 DO 154 1 = 1, NTT 1-034 

154 HH( n=HS( I ) 1-035 

GO TO 156 1-036 

153 DO 155 1=1, NTT 1-037 

155 HS( I )=HH( I ) 1-038 

WRITE OUTPUT TAPE 6,117 1-039 

WRITE OUTPUT TAPE 6 , 1 19 , { HH ( I ) ,XX ( I ) ,WX ( I ) , 1= 1 ,NTT ) 1-040 

156 CONTINUE 1-041 



101 F0RMAT(13A6) 




1-042 

102 F0RMAT(20H1 

ELASTIC SCATTERING/1H013A6 ) 

1-043 

103 F0RMAT(20H0 

INPUT TO PART 

3NE ) 


1-044 

104 FDRMAT(9H0 

AVI =1PE15.8* 

9H 

AWl =1PE15.8* 

1-045 

1 9H 

ARl =1PE15.8* 

9H 

AAl =1PE15.8) 

1-046 

105 F0RMAT<9H0 

DVl =1PE15.8* 

9H 

OWl =1PE15.8* 

1-047 

1 9H 

DRl =1PE15.8> 

9H 

DAI =1PE15.8) 

1-048 

144 F0RMAT(9H0 

NVl =115 * 

9H 

NWl =115 * 

1-049 

1 9H 

NRl =115 * 

9H 

NAl =115 ) 

1-050 

106 FDRMAT(9H0 

AV2 =1PE15*8* 

9H 

AW2 =1PE15.8* 

1-051 

1 9H 

AR2 =1PE15#8* 

9H 

AA2 =1PE15.8) 

1-052 

107 F0RMAT(9H0 

DV2 =1PE15.8* 

9H 

DW2 =1PE15.8* 

1-053 

1 9H 

DR2 =1PE15.8» 

9H 

DA2 =1PE15.8) 

1-054 

145 FDRMAT{9H0 

NV2 =115 * 

9H 

NW2 =115 * 

1-055 

1 9H 

NR2 =115 * 

9H 

NA2 =115 ) 

1-056 

108 F0RMAT(9H0 

AEF =1PE15.8» 

9H 

ARC =1PE15.8* 

1-057 

1 9H 

DEL =1PE15.8* 

9H 

ARR =1PE15.8) 

1-058 

109 F0RMAT(9H0 

AMP =1PE15.8* 

9H 

AZP =1PE15.8* 

1-059 

1 9H 

AMT =1PE15.8* 

9H 

AZT =1PE15.8) 

1-060 

110 F0RMAT{9H0 

ATI =1PE15.8* 

9H 

DTI =1PE15.8* 

1-061 

1 9H 

NTl =115 , 

9H 

NTT =115 ) 

1-062 

111 F0RMAT{9H0 

BBl =1PE15.8* 

9H 

BB2 =1PE15.8* 

1-063 

1 9H 

BB3 =1PE15«8* 

9H 

BB4 =1PE15.8 ) 

1-064 

112 F0RMAT(9H0 

BB5 =1PE15.8» 

9H 

BB6 =1PE15.8* 

1-065 

1 9H 

BBT =1PE15.8* 

9H 

BB8 =1PF15.8 ) 

1-066 

113 F0RMAT(9H0 

BB9 =1PE15.8* 

9H 

BIO =1PE15.8* 

1-067 

1 9H 

Bll =1PE15.8) 



1-068 

114 F0RMAT(9H0 

NNl =115 * 

9H 

NN2 =115 * 

1-069 

1 9H 

NN3 =115 * 

9H 

NN4 =115 ) 

1-070 

115 F0RMAT(9H0 

NN5 =115 , 

9H 

NN6 =115 * 

1-071 

1 9H 

NN7 =115 * 

9H 

NN8 =115 ) 

1-072 

160 FDRMAT(9H0INK0DF =115 , 

9H 

KDPLT =115 * 

1-073 

1 9H KDGRD =115 , 

9H 

KDCAL =115 ) 

1-074 

116 F0RMAT(8F10* 

0) 



1-075 

117 F0RMAT(16H1 

CM ANGLE (DEG) 

*20H 

CROSS SECTION (MB) ♦ 

1-076 

1 


21H 

ERROR WEIGHT FACTOR) 

1-077 

118 FDRMAT(8E10. 

0) 



1-078 

119 F0RMAT(0PF16 

• 3*1P2E20*8 ) 



1-079 

120 F0RMAT( 1515) 




1-080 

125 F0RMAT(6F10. 

0) 



1-081 

REDUCED MASS 




1-082 

X= i AMP*AMT ) / ( AMP+AMT ) 



1-083 

WAVE NUMBER 




1-084 

AEE = AEE*X/AMP 



1-085 

AKS = (2. 0*X*AEE)/41. 826134 


1-086 

AKK = SQRTF{AKS) 



1-087 

MESH SIZE 




1-088 

ADD=1000.0*DEL/AKK 



1-089 

ADD=0.001*INTF( ADD+0.5 ) 



1-090 

COULOMB PARAMETER 



1-091 

AZZ= (AZP*AZT*X*0. 034444017) /AKK 

1-092 

A2D=2.0*ADD 




1-093 

A3D=3.0*ADD 




1-094 

A6D=6*0*ADD 




1-095 

NORMALIZATION FACTOR FOR L 

=0 RADIAL WAVE FUNCTION 

1-096 

IF(AZZ) 121* 

122*121 



1-097 

122 ACC=1.0 




1-098 

GO TO 123 




1-099 

121 Y=6,2831853*AZZ 



1-100 

ACC=SQRTF( Y/ ( EXPF (Y)-l .0) ) 



1-101 

INTERVAL DOUBLING RADIUS 



1-102 
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123 IRR=ARR/A3D+0.50 1-103 

IRR=3*IRR 1-104 

ARR=ADD*FL0ATF( IRR) 1-105 

C MAXIMUM RADIUS 1-106 

Y=ARR+12.0*ADD 1-107 

X=AR1+AA1#BB1+BB2 1-108 

Z=AR2+AA2*BB1+BB2 1-109 

W=BB3+AZZ/AKK 1-110 

ARM=MAX1F(X,Y,Z*W) 1-111 

IRM=(ARM-ARR) /A6D+0.50 1-112 

IRM=3*IRM 1-113 

ARM=ARR+A2D*FL0ATF( IRM) 1-114 

IRM=IRR+IRM 1-115 

D12=ADD*ADD/12,0 1-116 

C MAXIMUM L VALUE 1-117 

LRM=AKK*ARM+BB4 1-118 

LRM=XMIN0F(LRM,50) 1-119 

DDS=D12*AKS/AEP 1-120 

X=AZZ*AZZ 1-121 

C NORMALIZATION FACTORS FOR RADIAL WAVE FUNCTION 1-122 

DO 124 I=1*LRM 1-123 

Y=I 1-124 

124 00(1 )=SORTF( Y*Y+X)/(Y*(2.0*Y+1.0) ) 1-125 

WRITE OUTPUT TAPE 6»131 1-126 

WRITE OUTPUT TAPE 6 * 132 •AKK • ADD ,AZZ , ACC 1-127 

WRITE OUTPUT TAPE 6 , 1 33 » ARR » ARM , I RR ♦ I RM 1-128 

WRITE OUTPUT TAPE 6»134»LRM 1-129 

131 F0RMAT(21H1 OUTPUT FRO*^ PROGRAM) 1-130 

132 F0RMAT(9H0 AKK =1PE15.8, 9H ADD =1PE15.8* 1-131 

1 9H AZZ =1PE15*8, 9H ACC =1PE15.8) 1-132 

133 F0RMAT(9H0 APR =1PE15.8* 9H ARM =1PE15.8» 1-133 

1 9H IRR =115, 9H IRM =115 ) 1-134 

134 F0RMAT(9H0 LRM =115 ) 1-135 

C PREPARATION FOR PROGRAM EIGHT 1-136 

II4=NN6 1-137 

K800-0 1-138 

RETURN 1-139 

END 1-140 
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SUBROUTINE TWO 

C PROGRAM TWO COULOMB PHASE SHIFTS AND WAVE FUNCTIONS 

COMMON AAl ,AA2 , ACC , ADD , AKK » AKS * ARC , ARM ,ARR »AR 1 , AR2 » ASS » ATI » AVI . 

1 AV2»AW1,AW2,AZD»AZZ»A2D»A3D»A6D> 

2 BB1,BB2,BB3,BB4,BB5*BB6*BB7,BB8*BB9,B10,B11»B12.BM,BP* 

3 CC,CN,CS* 

4 DAI ,DA2,DD*DDS,DR1 ,DR2,DT1 »DV1 ,DV2»DW1 ,DW2 ,D12* 

5 ERR»FF,GG*HH,HS, 

6 I 11*112,1 13,1 I4,ICC»INK0DE»IRM,IRR»I3I , 

7 KDCAL ,KDGRD ,KDPLT ,KODE2 ,K00E6 ,K0DE7 ,K0DE8 ,K800, 

8 LL,LLM,LRM, 

9 NA1X,NA2X,NN1 ,NN2 ,NN3 , NN4 , NN5 ,NN6 ,NN7 , NN8 ,NNN ,NR1X ,NR2X 
COMMON NTN,NTT,NT1 , NVl X , NV2X , NWl X , NW2X ,N3N ,PL , 

1 QQ,RR,SN,SS,TT,VC,VV,WW,WX,XC,XN,XX 

DIMENSION BM(51) ,BP(51) , 

1 CC(2,100) ,CN(90) ,CS(51) , 

2 DD(2,100) ,FF(2,51),G6(2,51) ,HH(90) ,HS(90) ,PL(51) , 

3 00(50) ,RR(500) ,SS(500) , SN ( 5 1) > TT ( 50 0 ) ,VC(500) ,VV(500) , 

4 WW(500) ,WX(90) , 

5 XC(90) ,XN(90) ,XX(90) 

DIMENSION AA( 100) ,AS( 51 ) ,3B(100) ,FA( 2,101 ) 

C COULOMB DHASE SHIFTS 

A=ATANF( AZZ/51.0) 

B=SQRTF(2601 .0+AZZ**2 ) 

C STIRLINGS FORMULA GIVE'^ AS (L + 1) FOR LARGE L 

AS(51)=50,5*A+AZZ*L0CF(B)-AZZ-SINF(A)/(12.0*B)+SINF(3. 

1 0*A) /{ 360*0*8**3 )-SINF(5.0*A)/(126C.0*B**5)+5INF(7.0*A)/{ 1680 

2 .0*B**7) 

DO 201 1=1,50 
J=51-I 

201 AS( J)=AS( J+1)-ATANF(AZZ/FL0ATF(J) ) 

ASS=2*0*AS( 1 )+AZZ*0,69314718 

DO 202 1=1,51 
A=2.0*AS( I ) 

CS( I )=CDSF(A) 

202 SN( I )=SINF(A) 

C COULOMB WAVE FUNCTIONS FOR L=0,1 BY ASYMPTOTIC SERIES 

203 DO 204 M=l,2 
DO 205 N=l,2 
L=M-l 

r=ARM+A2D*FL0ATF ( N-2 ) 

A=2.0*AZZ 
T = L 

B=T*( T+1.0 )+AZZ*AZZ 

C=2 .0*R*AKK 

U=1.0 

V=0.0 

X=1.0 

Y = 0.0 

DO 206 J=1,NN1 
D = J 

P=(A-AZZ/D)/C 
Q=(B/D+1.0-D)/C 
ZR=X*P - Y*0 
ZI=X*Q + Y*P 

DIV= ZR*ZR + ZI*ZI - X*X - Y*Y 
IF(DIV)230,207,207 
230 X=ZR 
Y=ZI 
U=U+X 

V = V+Y 
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n n 


Z* ( X**2+Y**2 ) / nJ**2+V**2 ) 2-042 

IF(2-BB5) 207 ♦206*206 2-043’ 

206 CONTINUE 2-044 

GO TO 209 2-045 

207 P-R*AKK-AZZ*L0GF(C)+AS(L+1) -1.5707963*1 2-046 

X=:COSF(P) 2-047 

Y«SINF(P) 2-048 

FF ( N *M ) =X*V+Y*U 2-049 

205 GG( N*M) =X*U-Y*V 2-050 

204 CONTINUE 2-051 

C CHECK ACCURACY OF RESULT WITH WRONSKIAN RELATION 2-052 

A=1.0/SORTF( 1.0+AZZ**2) 2-053 

B^FF ( 1 ♦ n *GG (1*2) -FF ( 1 ♦ 2 > *GG ( 1 ♦ 1 > -A 2-054 

CsFF ( 2 ♦ 1 ) *GG (2*2)-FF(2*2) *GG ( 2 ♦ 1) -A 2-055 

BaABSF(B) 2-056 

C:=ABSF(C) 2-057 

I F( B-l.E-5) 2081 *209*209 2-058 

2081 IF(C-l.E-5)208*209*209 2-059 

209 ARM=ARM+A6D*FL0ATF(NN2) 2-060 

IRM=IRM+3*NN2 2-061 

IF( I RM-500 >203*260,260 2-062 

260 WRITE OUTPUT TAPE 6*240 2-063 

240 F0RMAT(68H0 ARM EXCEEDS AVAILABLE STORAGE FOR POTENTIALS* PROCEED 2-064 

ITO NEXT CASE) 2-065 

IF(KDPLT)261*262*261 2-066 

261 READ INPUT TAPE 7*241*PL0T 2-067 

241 FORMAT { FI 0.0 ) 2-068 

262 K0DE2=1 2-069 

RETURN 2-070 

COULOMB WAVE FUNCTIONS FOR L MORE THAN 1 BY RECURRENCE RELATIONS 2-071 

GET COEFFICIENTS FOR THE RECURRENCE RELATIONS 2-072 

208 L=LRM+10 2-073 

DO 210 I=1*L 2-074 

BB( I )=AZZ/FLOATF( I ) 2-075 

210 AA( I ) = SQRTF( 1.0+BB( n**2) 2-076 

J=L-1 2-077 

DO 211 N=l*2 2-078 

R=AKK*(ARM+A2D*FL0ATF(N-2 ) ) 2-079 

DO 212 I=1*J 2-080 

212 CC(N*I)=FL0ATF(2*I+1)/R+BB(I)+BB(I+1) 2-081 

211 CONTINUE 2-082 

C RECURR UP TO GET THE IRREGULAR COULOMB WAVE FUNCTIONS 2-083 

J=LRM+1 2-084 

DO 213 N=l*2 2-085 

DO 214 I=3*J 2-086 

214 GG( N* I) = (CC{ N* 1-2 )*GG ( N * I-l )-AA ( 1-2 ) *GG { N ♦ 1-2 ) ) /AA( I-l ) 2-087 

213 CONTINUE 2-088 

C RECURR DOWN TO GET THE REGULAR COULOMB WAVE FUNCTIONS 2-089 

215 J=L+1 2-090 

FAC 1*J)=:0.0 2-091 

FA(2*J)=0.0 2-092 

FAC 1,L)=0.1 2-093 

FAC2*L)=0.1 2-094 

K=L-1 2-095 

DO 216 N=l*2 2-096 

DO 217 I=1*K 2-097 

J=L-I 2-098 

FACN*J)=CCCCN*J)*FACN*J+1 )-AAC J+l)*FACN*J+2) ) /AAC J) 2-099 

IF( FA CN*J)-1.E30) 217,231,231 2-100 

231 IFC J-LRM+1 >232*233,233 2-101 

233 FACN*J)=FACN*J)*l.E-30 2-102 
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FA(N,J+l)=FA(N,J+l)*l.E-30 2-103 

GO TD 217 2-104 

232 JK=LRM+1 2-105 

DO 234 IJ=J,JK 2-106 

234 FA(N,IJ)=FA{N*IJ)*l.E-30 2-107 

217 CONTINUE 2-108 

216 CONTINUE 2-109 

C RENORMALIZE THE REGULAR COULOMB WAVE FUNCTIONS 2-110 

C=FF(1,1)/FA(1,1) 2-111 

D=FF{2tl)/FA(2,l) 2-112 

DO 218 I=1*L 2-113 

FA( 1 ,I )=C*FA(1,I ) 2-114 

218 FA(2»n=D*FA(2*I ) 2-115 

C CHECK ALL COULOMB WAVE FUNCTIONS WITH WRONSKIAN RELATION 2-116 

M=LRM+1 2-117 

DO 219 I=3»M 2-118 

D = FA(1 * I-1)*GG(1 * I )-FA( 1 *I )*GG( 1 » I -1 ) -1 . 0/ AA ( I-l) 2-119 

IF( ABSF(D)-1.E-5)219»219»221 2-120 

219 CONTINUE 2-121 

GO TO 220 2-122 

C RECURRENCE RELATION HAS FAILED INCREASE L AND TRY AGAIN 2-123 

221 J=L+1 2-124 

L=L+10 2-125 

IF(L-100)222»222*223 2-126 

223 WRITE OUTPUT TAPE 6,229 2-127 

229 FDRMAT(74H0RECURRENCE RELATIONS FDR COULOMB WAVE FUNCTION HAVE FAI 2-128 

ILED, INCREASE RMAX/ ) 2-129 

GO TO 209 2-130 

222 DO 224 I=J,L 2-131 

C ADDITIONAL COEFFICIENTS FDR THE RECURSION RELATIONS 2-132 

BB{ I )=AZZ/FLDATF( I ) 2-133 

224 AA( I )=SQRTF( 1.0+BB( I )** 2 ) 2-134 

M=J-1 2-135 

K=L-1 2-136 

LLM=K 2-137 

DO 225 N=l,2 2-138 

r=AKK*{ARM+A2D*FLDATF(N-2 )) 2-139 

DO 226 I=M,K 2-140 

226 CC (N, n =FLDATF ( 2*1 + 1 ) /R + BB( I ) +BB ( I + l ) 2-141 

225 CONTINUE 2-142 

GO TO 215 2-143 

220 DO 228 N=l,2 2-144 

DO 227 1^1,51 2-145 

227 FF(N,I )=FA(N,I ) 2-146 

228 CONTINUE 2-147 

250 F0RMAT(15H0 FINAL VALUE 5X,5HARM = 1 PE 1 5 . 8 , 5 X , 5H I RM =14) 2-148 

WRITE OUTPUT TAPE 6 , 250 , ARM , I RM 2-149 

K0DE2=2 2-150 

RETURN 2-151 

END 2-152 
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SUBROUTINE THREE 3-000 

C PROGRAM THREE COULOMB POTENTIAL AND R SQUARED 3-001 

COMMON AA1,AA2,ACC,ADD,AKK»AKS,ARC,ARM»ARR,AR1,AR2,ASS*AT1>AV1» C-Ol 

1 AV2 ,AW1,AW2,AZD,AZZ,A2D,A3D,A6D, C~02 

2 BB1,BB2,BB3,BB4,BB5,BB6,BB7,BB8,BB9,B10*B11,B12,BM*BP» C-03 

3 CC*CN»CS, C-04 

4 DAI ,DA2,DD»DDS*DR1*DR2,DT1,DV1,DV2»DW1 ,DW2*D12, C-05 

5 ERR*FF,GG»HH,HS» C-06 

6 II1,II2,II3,II4,ICC»INK0DE>IRM,IRR,I3I» C-07 

7 KDCAL,KDGRD,KDPLT,KODE2»KODE6*KODE7*KODE8»K800» C-08 

8 LL,LLM,LRM, C-09 

9 NA1X,NA2X»NN1 »NN2 »NN3 , NN4 » NN5 * NN6 f NN7 , NN8 »NNN ,NR IX *NR2X C-10 

COMMON NTN,NTT,NTl,NVlX*NV2X,NWlX»NW2X»N3NfPL* C-11 

1 QQ»RR,SN,SS,TT,VC»VV»WW,WX*XC.XN*XX C-12 

DIMENSION BM( 51) *BP(51 ) » D-01 

1 CC(2,10C) ,CN(90) ,CS(51) * D-02 

2 DD(2»100),FF(2,51)*GG(2»51)»HH(90),HS(90),PL(51)> D-03 

3 QQ( 50) »RR ( 500) »SS ( 500) f 5N( 51 ) »TT( 500) *VC( 500) »VV( 500 ) * D-04 

4 WW(500) »WX(90) » D-05 

5 XC(90) *XN(90 ) ♦XX(90) D-06 

C PUT RC ON THE MESH 3-002 

R = ARC - ARR 3-003 

IF(R) 301*301»302 3-004 

301 I=ARC/A3D+0.5 3-005 

J=3*I 3-006 

ARC=ADD*FLOATF ( J ) 3-007 

GO TO 303 3-008 

302 I=R/A6D+0.5 3-009 

J=3*I 3-010 

ARC=ARR+A2D*FLDATF ( J) 3-011 

J=IRR+J 3-012 

303 WRITE OUTPUT TAPE 6.304*ARC,J 3-013 

304 F0RMAT(15H0 FINAL VALUE5X ,5HARC = IPE 1 5 . 8 , 5X » 5H I RC =14//) 3-014 

ICC=J 3-015 

C CALCULATE THE COULOMB POTENTIAL 3-016 

AZD=2.0*AKK*AZZ 3-017 

W=D12*AZD 3-018 

X=D12*AKS 3-019 

Y=W*1. 5/ARC 3-020 

Z = Y*0. 3 3333333 /{ARC**2 ) 3-021 

Y=X-Y 3-022 

R=0.0 3-023 

K=XMIN0F( J,IRR) 3-024 

DO 305 1=1, K 3-025 

R=R+ADD 3-026 

RR(I)=R*R 3-027 

305 VC( n=Y+Z*RR(I ) 3-028 

IF(J-IRR) 306,307,308 3-029 

306 K=J+1 3-030 

DO 309 I=K,IRR 3-031 

R=R+ADD 3-032 

RR{I)=R*R 3-033 

309 VC( I )=X-W/R ‘ 3-034 

K=IRR+1 3-035 

X=4.0*X 3-036 

W=4.0*W 3-037 

GO TO 310 3-038 

308 K=IRR+1 3-039 

Y=4.0*Y 3-040 

Z=4.0*Z 3-041 

DO 311 I=K,J 3-042 

R=R+A2D 3-043 

RR( I )=R*R 3-044 

311 VC( I )=Y+Z*RR( I) 3-045 

307 K=J+1 3-046 

W=4.0*W 3-047 

X=4,0*X 3-048 

310 DO 312 I=K,IRM 3-049 

R=R+A2D 3-050 

RR( I)=R*R 3-051 

312 VC( I )=X-W/R 3-052 

RETURN 3-053 

END 3-054 
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SUBROUTINE FOUR 4-000 

C PROGRAM FOUR OPTICAL POTENTIAL FORM FACTORS AND POWER SERIES 4-001 

C COEFFICIENTS 4-002 

COMMON AA1,AA2,ACC,ADD»AKK,AKS,ARC»ARM,ARR ,AR 1 , AR2 » ASS *AT1 •AVI* C-01 

1 AV2,AW1 ,AW2,AZD*AZZ*A2D*A3D*A6D* C-02 

2 BB1*BB2 *BB3,BB4*BB5*BB6*BB7*BB8 *BB9 *B 1 0 * B1 1 * B 1 2 * BM • BP ♦ C-03 

3 CC*CN*CS* C-04 

4 DAI *DA2,DD*DDS*DR1 *DR2*DT1 *DV1 *DV2*0W1 *DW2 *D12* C-05 

5 ERR>FF*GG»HH*HS* C-06 

6 I I1,I 12*1 13* 1 14* ICC*INKODE*IRM, IRR, I3I * C-07 

7 KDCAL*KDGRD*KDPLT*KODE2*KODE6*KODE7*KODE8*K800* C-08 

8 LL*LLM*LRM» C-09 

9 NA1X*NA2X»NN1*NN2*NN3*NN4*NN5*NN6*NN7*NN8*NNN*NR1X*NR2X C-10 

COMMON NTN*NTT*NT1*NV1X*NV2X*NW1X*NW2X*N3N*PL * C-11 

1 QQ*RR*SN*SS*TT*VC*VV*WW*WX*XC*XN*XX C-12 

DIMENSION BM(5l)*BP(51)* D-01 

1 CC(2*100) ,CN(90) *CS(51) * D-02 

2 D0(2*100) *FF(2*51 ) ,GG(2*51) *HH(90) *HS(90) *PL(51) * D-03 

3 0Q( 50) *RR( 500) *SS(500) *SN( 51) *TT(500) *VC( 500) *VV{ 500) * D-04 

4 WW( 500) *WX ( 90) * D-05 

5 XC(90) *XN(90) >XX(90) D-06 

DIMENSION AC(IOO) D-07 

C CALCULATION OF SAXON FORM FACTOR 4-003 

401 A=EXPF( ADD/AAl ) 4-004 

B=EXPF(-AR1/AA1 ) 4-005 

DO 402 1=1. IRR 4-006 

B=B*A 4-007 

402 SSI I ) = 1.0/(1,0+B) 4-008 

A=A*A 4-009 

J=IRR+1 4-010 

DO 403 I=J*IRM 4-011 

B=B*A 4-012 

SSI I ) = 1.0/tl.0+B) 4-013 

IFI SSI I )-BB6 ) 404*403.403 4-014 

403 CONTINUE 4-015 

GO TO 405 4-016 

404 NNN=XMAXOFI I .lRR+6) 4-017 

J=I+1 4-018 

DO 406 I=J.IRM 4-019 

406 SSI I )=0,0 4-020 

GO TO 407 4-021 

405 NNN=IRM+1 4-022 

407 IFINN3) 408*409*408 4-023 

C CALCULATION OF SURFACE INTERACTION FORM FACTOR 4-024 

408 X=ADD/AA2 4-025 

Y=-AR2/AA2 4-026 

Z=L0GFIBB6) 4-027 

DO 410 I=1*IRR 4-028 

Y=Y+X 4-029 

IFIY-Z) 410*412*412 4-030 

410 TTI I )=0.0 4-031 

GO TO 413 4-032 

412 NTN=I 4-033 

Y=EXPF(Y-X) 4-034 

X=EXPFtX) 4-035 

DO 414 J=I *IRR 4-036 

Y=Y*X 4-037 

414 TTIJ)=4.0/I2.0+Y+1*0/Y) 4-038 

X=X*X 4-039 

J=IRR+1 4-040 

415 DO 416 I=J*IRM 4-041 
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n n 


n = i 

Y=Y*X 

TT( I ) = ^.0/(2«0 + Y+1,0/Y) 

IF(TT( I )-BB6) 417*A16»416 

416 CONTINUE 
I = I I 

417 NNN=XMAXOF(NNN»I ) 

IF( I-IRM) 418»409»418 

418 K=I+1 

DO 420 J=KfIRM 
420 TT(J)=0.0 
GO TO 409 
413 X=2.0*X 
J=IRR+1 

DO 422 I=J,IRM 
Y = Y + X 

IF( Y-Z ) 422,423*423 

422 TT( I ) = 0.0 
GO TO 409 

423 NTN=I 
Y=EXPF( Y-X) 

X=EXPF(X) 

J=I 

GO TO 415 
409 A=NNN-IRR 
N3N=A/6.+,5 
NNN=IRR+3*N3N 
I3I=( IRM-NNN) /3 

CALCULATION OF POWER SERIES COEFFICIENTS 
FIRST CALCULATE THE SAXON WELL COEFFICIENTS 
A=EXPF(-AR1/AA1) 

CC(1,1)=1.0/(1.0+A) 

B=-A*CC(1*1) 

AC( 1) = 1.0 * 

CC( 1 ,2 )=-A*CCt 1 ,1 )**2/AAl 
AC(2)=1.0 

D=CC( 1*1 )/(2.0*AA1*AA1 ) 

CC( 1,3)=CC(1,1)*CC(1,2)*( 1.0-A)/(2.0*AA1) 

DO 451 I=3*NN4 
AC( I ) = 1.0 
DO 452 J=3,I 
K=I+2-J 

452 ACC K)=AC(K)#FL0ATF{K)+AC{ K-1) 

K=I + 1 

D = D/{ AA1*FLDATF( I ) ) 

IF( D-1 .E“30) 463*463*464 
464 E=l. 

F = 0.0 

DO 453 J=1*I 
E=E*B*FLOATFC J) 

IFC ABSFC E)-l.E-30) 461*461*453 
451 E=0. 

453 F=F+AC(J)*E 
451 CCC 1,K)=F*D 
463 DO 467 I I=K*LLM 
467 CCC 1*1 I )=0. 

459 CONTINUE 
RETURN 
END 
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SUBROUTINE FIVE 

C PROGRAM FIVE THE EFFECTIVE KINETIC ENERGY AND COMBINED 

C POWER SERIES COEFFICIENTS 

COMMON AAl *AA2,ACC»ADD»AKK*AKS*ARC *ARM ,ARR * AR 1 * AR 2 * ASS > ATI » AVI » 

1 AV2 »AW1 ,AW2 ,AZD*AZZ *A2D»A3D*A6D » 

2 BBl *BB2 ,BB3,BBA,BB5>BB6 »BB7 *BB8 *BB9 ,BlO*Bll»B12*BM*BP* 

3 CC»CN*CS> 

A DAI *DA2 ,DD*DDS,DR1 *DR2 *DT1 *DV1 *DV2 »DW1 *DW2 *D12» 

5 ERR,FF*GG»HH,HS» 

6 I I1»I 12,1 I3,IIA,ICC,INK0DE.IRM.IRR»I3I , 

7 KDCAL*KDGRD,KDPLT,KQDE2»KDDE6*KODE7,KODE8»K800, 

8 LL*LLM»LRM» 

9 NA1X*NA2X*NN1 ,NN2 ,NN3 , NN4 ,NN5 ,NN6 ,NN7 , NN8 »NNN ,NR IX * NR2X 
COMMON NTN,NTT,NT1,NV1X»NV2X.NW1X*NW2X »N3N»PL* 

1 QO,RR,SN»SS,TT,VC,VVtWW*WX»XC*XN*XX 

DIMENSION BM( 51 ) *BP{ 51 ) * 

1 CC<2,100) ,CN{90) »C5(51) » 

2 DD( 2 ,100) ,FF( 2 ,51 ) »GG( 2 »51 ) *HH(90) *HS (90) >PL( 51 ) * 

3 QQ( 50) ,RR{ 500) ,SS( 500) , SN ( 5 1 ) , TT ( 50 0 ) ,VC( 500) ,VV( 500) , 

4 WW(500) ,WX(90) , 

5 XC(90) ,XN(90) ,XX(90) 

C CONTRIBUTION TO E.K,E« DUE TO K.E., SAXON POTENTIAL, AND 

C COULOMB POTENTIAL 

500 V=-AV1*DDS 
W=-AW1*DDS 

DO 501 1 = 1, IRR 
VV( I )=VC( I )-V*SS { I ) 

501 WW( I )=-W*SS( I ) 

J=IRR+1 

V=4«0*V 

W=4.0-«-W 

DO 502 I=J,NNN 
VV( I )=VC( I )-V*SS( I ) 

502 WW( I )=-W*SS( I ) 

J=NNN+1 

DO 503 I=J,IRM 
WW( I )=0*0 

503 VV( I )=VC( I ) 

IF( NN3 ) 505,504,505 

C CONTRIBUTION TO E.K.E. DUE TO SURFACE INTERACTION POTENTIAL 

505 V=-AV2*DDS 
W=-AW2*DDS 

IF(NTN-IRR) 507,507,508 

507 DO 506 I=NTN,IRR 
VV( I )=VV( I )-V*TT( I ) 

506 WW( I )=WW(I )-W*TT ( I ) 

GO TO 516 

508 J=NTN 

GO TO 515 
516 J=IRR+1 
515 V=4*0*V 
W=4.0*W 

DO 509 I=J,NNN 
VV( I )=VV( I )-V*TT( I ) 

509 WW( I )=WW( I )-W*TT ( I ) 

C CALCULATE COMBINED POWER SERIES COEFFICIENTS 

504 A=-DDS/D12 
S=-A*AV1 
T=-A*AW1 

511 DO 514 1=1, NN4 

DD(1,I)=S*CC(1,I) 

514 DD(2,I )=T*CC( 1,1 ) 

513 CONTINUE 
RETURN 
END 
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SUBROUTINE SIX 6-000 

C PROGRAM SIX THE RADIAL WAVE FUNCTIONS AND ASSOCIATED SCATTERING 6-001 

C amplitudes 6-002 

COMMON AA1*AA2,ACC,ADD,AKK»AK5»ARC,ARM,ARR,AR1,AR2»ASS*AT1*AV1* C-Ol 

1 AV2,AW1,AW2,AZD»AZZ,A2D,A3D»A6D» C-02 

2 BB1»BB2,BB3,BBA»BB5»BB6»BB7,BB8»BB9»B10»B11*B12,BM*BP» C-03 

3 CC*CN,C5, C“OA 

4 DA1»DA2,DD»DDS»DR1 ,DR2»DT1 *DV1 ,DV2 »DWI ,DW2 ,D12 * C“05 

5 ERR,FF,GG»HH,HS» C-06 

6 IIl,II2,II3,Il4,lCC,INKaDEf IRM,IRR,I3I , C-07 

7 KDCAL,KDGRD,KDPLT,KODE2,KODE6,KODE7,KODE8*K800> C-08 

8 LL,LLM,LRM, C-09 

9 NA1X»NA2X»NN1»NN2,NN3*NN4>NN5*NN6»NN7»NN8»NNN*NR1X»NR2X C-10 

COMMON NTN,NTT,NT1 »NV1X >NV2X*NW1X *NW2X »N3N*PL » C-11 

1 QQ»RR,SN*SS,TT,VC,VV,WW»WX*XCfXN»XX C-12 

DIMENSION BM(5l)»BP(51)t D-01 

1 CC(2,100) ,CN(90) ,CS(51) » D-02 

2 DD(2,100) ,FF(2,51) »GG(2*51) ,HH(90) »HS(90) »PL(51) » D-03 

3 QQ( 50) »RR ( 500) ,SS( 500) »SN( 51 ) »TT(500) *VCt 500) *VV( 500) * D-04 

4 WW( 500) »WX{90) , D-05 

5 XC(90) »XN(90) >XX(90) D-06 

DIMENSION BB(2,100) D-07 

600 P=0*0 6-003 

MMM=0 6-004 

ICC=ICC 6-005 

C GET STARTING VALUE FOR I 6-006 

LL=0 6-007 

601 Q={2.*P+1.) 6-008 

I3J=I3I 6-009 

R=B11*P*ADD 6-010 

615 RT=P*ADD/3,5 6-011 

IF(R-RT)614t651*651 6-012 

651 S=( ABSF(R-ARR) + (R-ARR))/2.0 6-013 

IF(S) 602*602*603 6-014 

602 I=R/A3D+0.5 6-015 

I=XMAXOF{ 1,1) 6-016 

1=3*1 6-017 

R=ADD*FLOATF( I ) 6-018 

GO TO 604 6-019 

603 I=S/A6D+0.5 6-020 

1=3*1 6-021 

R=ARR+A2D*FLDATF( I ) 6-022 

I=IRR+I 6-023 

C INITIAL MESH SIZE 6-024 

604 IF(I-IRR) 605*606*606 6-025 

605 E=ADD 6-026 

L=1 6-027 

GO TO 607 6-028 

606 E=A2D 6-029 

L=2 6-030 

607 IF(R-ARC) 608*608,609 6-031 

608 X=0,0 6-032 

Y=1.0 6-033 

GO TO 610 6-034 

609 X=1 .0 6-035 

Y=0.0 6-036 

C GET FIRST COUPLE TERMS OF POWER SERIES 6-037 

610 B=1,5*AZD*Y/ARC 6-038 

C=X*AZD 6-039 

A=B-AKS 6-040 

S=R-E 6-041 
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B=0.33333333*B/RR( ICC) 6-042 

BB(1,1)=1,0 6-043 

BB(2*1)=0.0 6-044 

BB( 1,2)=C/(Q+1.0) 6-045 

BB(2.2)=0.0 6-046 

T=2.0*(0+2.0) 6-047 

BB ( 1 , 3 ) = ( A-DO ( 1 • 1 ) +C*BB ( 1 * 2 ) ) /T 6-048 

F=R*R 6-049 

G=S*S 6-050 

Z=F*R 6-051 

W=G*S 6-052 

BB( 2*3)=(-DD(2»l )/T) 6-053 

T=3.0*(Q+3*0) 6-054 

BBC 1*4)=(-DD( 1*2)+C*BB( 1*3)+A*BB( 1 » 2 ) -DD ( 1 • 1 ) *BB ( 1 »2 ) ) /T 6-055 

BB(2»4)=(-DD{2*2 )+C*BB(2»3)-DD(2*l )*BB(lf2) )/T 6-056 

X = BB ( 1 » 1 ) +BB ( 1 » 3 ) *F+BB ( 1 » 4 ) *Z+BB ( 1 > 2 ) »R 6-057 

U= BB ( 1 , 1 ) +BB (1,3) *G+BB (1,4) *W+BB ( 1 , 2 ) *S 6-058 

Y = BB( 2 ,3 ) *F+BB ( 2 *4)*Z 6-059 

650 V=BB(2,3)*S+BB(2»4)*W 6-060 

C GET ADDITIONAL TERMS TO POWER SERIES BY RECURSION 6-061 

DO 611 N=5,NN4 6-062 

W=W*S 6-063 

Z-Z*R 6-064 

T=N-1 6-065 

T=T*(T+Q) 6-066 

G=0.0 6-067 

K=N-2 6-068 

H=0,0 6-069 

DO 612 J=1,K 6-070 

M=N-J-1 6-071 

G=G+DD( 1,J)*BB( 1,M)-DD(2 ,J)*BB(2,M) 6-072 

612 H=H+DD( 1,J)*BB(2*M)+DD(2,J)*BB(1 ,M) 6-073 

BB ( 1 , N ) = ( A*BB ( 1 , N-2 ) -B*BB ( 1 ,N-4 ) -G+C*BB ( 1 , N-1 ) ) /T 6-074 

BB { 2 , N ) = ( A*BB ( 2 ,N-2 ) -B*BB ( 2 ,N-4 ) -H+C*BB ( 2 ,N-1 ) ) /T 6-075 

G=BB(1,N)*Z 6-076 

H=BB(2,N)*Z 6-077 

X=X+G 6-078 

Y=Y+H 6-079 

U=U+BB(1,N)*W 6-080 

V=V+BB(2,N)*W 6-081 

IF(XM0DF(N,2) ) 639, 640, 639 6-082 

639 IF(MN8-N)613,613,640 6-083 

640 F=(G**2+H**2) /(X**2+Y**2) 6-084 

IF(F-BB7) 613,611*611 6-085 

611 CONTINUE 6-086 

R=R-6.0*E 6-087 

MMM=MMM+1 6-088 

IF(R-A3D) 614,615,615 6-089 

614 WRITE OUTPUT TAPE 6,641,MMM,P 6-090 

641 F0RMAT(39H0P0WER SERIES FAILED TO CONVERGE IN S I X/5H0MMM= I 5 , 5X • 2HP 6-091 

1=F10.6) 6-092 

IF ( KDPLT) 661,662,661 6-093 

661 READ INPUT TAPE 7, 663, PLOT 6-094 

663 FORMAT(FIO.O) 6-095 

662 K0DE6=1 6-096 

RETURN 6-097 

C NORMALIZE WAVE FUNCTIONS 6-098 

613 A=AKK*R 6-099 

B=AKK*S 6-100 

F=A*ACC 6-101 

6=B*ACC 6-102 
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IF(LU 616,616*617 

617 DD 618 K=1,LL 
F=F*A*00( K) 

618 G=G*B*QQ(K) 

616 U=U*G 

V = V*G 
X = X*F 
Y=Y*F 

NOW integrate the diff equation 

SET UP INITIAL VALUES 
H=D12*P*(P+1*0) 

F=1.0+VV( I)-H/RR( I ) 

G=WW( I ) 

JJ=I-L 

C=1.0+VV( JJ)-H/RR(JJ) 

D=WW( JJ) 

ASSIGN 622 TO II 
IF(L-l) 619,619,620 

619 J=( IRR-I ) /3 

C INTEGRATION WITH A COMPLEX POTENTIAL 

624 DD 621 K=1,J 
1 = 1+1 

A=1.0+VV( I )-H/RR{ I ) 

B=WW( I ) 

Z=A*A+B*B 

0=12. 0-10. 0*F 

XP= { 0*X-C*U+10. 0*G*Y+D*V ) 

XM= (Q*Y-C*V-10.0*G*X-D*U) 
S=(A*XP+B*XM)/Z 
T=t A*XM-B*XP)/Z 
1 = 1 + 1 

C=1.0+VV( I )-H/RR( I ) 

D=WW( I ) 

Z=C*C+D*D 

Q=12.0-10.0*A 

XP= (Q*S-F*X+10.0*B*T+G*Y ) 

XM= (Q*T-F*Y-10.0*B*S-G*X ) 
U=(C*XP+D*XM) /Z 
V=(C*XM-D*XP)/Z 
1 = 1+1 

F=1.0+VV( I )-H/RR( I ) 

G=WW( I ) 

Z=F*F+G*G 

Q=12.0-10.0*C 

XP= (0*U-A*S+10.0*D*V+B*T ) 

XM= {Q*V-A*T-10.0*D*U-B*S ) 
X=(F*XP+G*XM)/Z 

621 Y=(F*XM“G*XP)/Z 

GO TO 11,(622,623) 

C CHANGE OF INTERVAL SIZE 

622 U=S 

J = N3N 

V = T 

F=4.0*F“3.0 
C=4,0*A~3.0 
G=4.0*G 
D=4.0*B 
GO TO 626 

620 J=(NNN-I )/3 

626 ASSIGN 623 TO II 
H=4.0*H 
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IF(J)63 1*623, '624 6-164 

631 I3J=I3I+J 6-165 

IF( I3J) 627, 632 *623 6-166 

C INTEGRATION WITH A REAL POTENTIAL 6-167 

623 DO 625 K=1,I3J 6-168 

1=1+1 6-169 

A = 1.0+W( I )-H/RR( I ) 6-170 

0=12. 0-10. 0*F 6-171 

S=(Q*X-C»U)/A 6-172 

T=(Q*Y-C*V)/A 6-173 

1=1+1 6-174 

C=1.0+VV( I )-H/RR( I ) 6-175 

0=12. 0-10. 0*A 6-176 

U=(Q*S-F*X)/C 6-177 

V=(0*T-F*Y) /C 6-178 

1=1+1 6-179 

F=1 .0+VV( I )-H/RR( I ) 6-180 

0=12. 0-10. 0*C 6-181 

X={0*u-A*S)/F 6-182 

625 Y=(0*V-A*T)/F 6-183 

C CALCULATION OF THE SCATTERING AMPLITUDES 6-184 

632 LL=LL+1 6-185 

A=X*FF( 1 ,LL)-U*FF (2,LL) 6-186 

B=Y*FF(1,LL)-V*FF(2,LL) 6-187 

C=X*GG( 1,LL)-U*GG(2,LL ) 6-188 

D=Y*GG( 1 ,LL) -V*GG ( 2 ,LL ) 6-189 

E= { A+D )**2+( B-C) **2 6-190 

BP(LL)=( ( (D+A)*( D-A)+(C+BI*<C-B) )/E)-1.0 6-191 

BM{ LL) =(-2.0)*(A*C+B*D) /E 6-192 

C SHALL WE GO AROUND AGAIN 6-193 

IF( ABSF(BP(LL) ) + ABSF(BM(LD) - BB8 ) 627,628,628 6-194 

628 IF(LL-LRM) 629,629,627 6-195 

629 P=P+1.0 6-196 

GO TO 601 6-197 

630 F0RMAT(9H0 MMM =115, 9H LRM =115, 9H LMAX =5X,2PE10.3) 6-198 

627 WRITE OUTPUT TAPE 6 , 630 ,MMM , LRM , P 6-199 

K0DE6=2 6-200 

RETURN 6-201 

END 6-202 
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SUBROUTINE SEVEN 7-000 

C PROGRAM SEVEN THE ELASTIC CROSS SECTION 7-001 

COMMON AAl *AA2 , ACC , ADD > AKK , AKS » ARC ♦ ARM *ARR *AR1 »AR2 » ASS » ATI »AV1 1 C-Ol 

1 AV2»AW1 *AW2,AZD*AZZ,A2D*A3D»A6D* C-02 

2 BB1*BB2,BB3*BBA*BB5*BB6»BB7»BB8*BB9»B10*B11»B12*BM,BP» C~03 

3 CC»CN*CS* C-04 

4 DAI »DA2fDD»DDS>DRl »DR2 »DT1 »DV1 *DV2»DW1 »DW2 »D12» C-05 

5 ERR*FF*GG*HH*HS» C-06 

6 II1,II2,I I3,II4,ICC,INK0DE,IRM*IRR,I3I * C-07 

7 KDCAL »KDGRD,KDPLT ,K0DE2 .K0DE6 »1CDDE7 ,K0DE8 »K800 * C“08 

8 LL,LLM,LRM» C-09 

9 NAlX*NA2X*NNl*NN2*NN3>NN4»NM5tNN6»NN7.NN8»NNN,NRlX»NR2X C-10 

COMMON NTN,NTT,NT1,NV1X>NV2X»NW1X»NW2X»N3N,PL • C-11 

1 QQ,RR,SN.SS,TT,VC»VV.WW*WX*XC*XN»XX C-12 

DIMENSION BM(51)*BP(51)» D“01 

1 CC(2»100) ,CN(90)»CS(51) * D-02 

2 DD(2»100) ,FF(2*51 ) >GG( 2»51 ) *HH(90) »HS(90) *PL( 51 ) * D“03 

3 QQ( 50) >RR(500) ,SS( 500) » SN ( 5 1 ) * TT ( 500 ) .VC ( 500) .VV( 500) . D“04 

4 WW(500) >WX(90) . D-05 

5 XC(90) .XN(90) .XX(90) D-06 

736 IF(NN5) 708.708.732 7-002 

732 K0DE7=2 7-003 

719 ASSIGN 702 TO KK 7-004 

ASSIGN 717 TO JJ 7-005 

ERR=0.0 7-006 

N=NTT 7-007 

706 DO 701 J=1,N 7-008 

X=COSF(HH(J) *0.017453293 ) 7-009 

A=(1.0-X) 7-010 

B=AZZ/A 7-011 

C=A5S-AZZ*L0GF(A) 7-012 

E=-B*COSF(C) 7-013 

F=-B*SINF(C) 7-014 

PL(1)=1.0 7-015 

PL(2)=X 7-016 

G=0.5*(BP(1)*SN(1 )+BM( 1)*CS(1 ) ) +1 . 5*X* ( BP ( 2 ) *SN ( 2 ) +BM ( 2 ) *CS ( 2 ) ) 7-017 

H= 0.5*(BM(1 )*SN( 1)-BP( 1)*CS( 1) ) + 1 . 5*X* ( BM ( 2 ) *SN ( 2 ) -BP ( 2 ) *CS ( 2 ) ) 7-018 

DO 718 K=3.LL 7-019 

D=K-1 7-020 

PL( K)= ( (2.0*D-1 .0)*X*PL( K-1 )-(D-1.0)*PL(K-2 ) ) /D 7-021 

A=PL(K)*(D+0.5) 7-022 

G = G+A*(BP (K)*SN(K)+BM(K.)*CS(K) ) 7-023 

718 H=H+A*(BM(K)*SN(K)-BP{K)*CS(K) ) 7-024 

XN< J) = ( ( E+G)**2+( F+H )**2 ) *10.0/AKS 7-025 

GO TO KK. (702.713) 7-026 

713 XC ( J ) = ( E**2+F**2 ) *10 . 0/AKS 7-027 

IF(XC( J) ) 721 .720.721 7-028 

720 CN(J)=0. 7-029 

GO TO 722 7-030 

721 CN( J)=XN( J)/XC( J) 7-031 

722 CONTINUE 7-032 

WRITE OUTPUT TAPE 6 .71 1 .HH( J) *XN( J ) .XC ( J ) .CN ( J ) 7-033 

GO TO 701 7-034 

702 ERR=ERR+WX( J)*( (XN( J)-XX( J) )/XX( J) )**2 7-035 

701 CONTINUE 7-036 

GO TO JJ. (717,715) 7-037 

703 FDRMAT(9H AVI =1PE15.8. 9H AWl =1PE15.8. 7-038 

1 9H ARl =1PE15.8, 9H AAl =1PE15.8) 7-039 

704 F0RMAT(9H AV2 =1PE15.8. 9H AW2 =1PE15*8. 7-040 

1 9H AR2 =1PE15.8. 9H AA2 =1PE15.8) 7-041 

705 F0RMAT(9H ERROR =1PE15.8) 7-042 
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717 WRITE OUTPUT TAPE 6 * 70 3 , AVI » AWl , AR 1 , AA 1 7-0A3 

WRITE OUTPUT TAPE 6 , 704 , AV2 , AW2 » AR 2 , AA2 7-044 

WRITE OUTPUT TAPE 6, 705, ERR 7-045 

RETURN 7-046 

710 F0RMAT(45H1 CALCULATED VALUES OF ELASTIC CROSS SECTIONS 7-047 

1 /16H0 CM ANGLE (DEG),18H EL XSCTN (MB) 7-048 

2 18H COUL XSCTN (MB) ,18H EL/CDUL 7-049 

711 F0RMAT(F16.3,1PE18.8,1PE18.8,1PE18.8) 7-050 

708 IF(KDGRD) 731,730,731 7-051 

731 K0DE7=3 7-052 

GO TO 719 7-053 

730 ASSIGN 713 TO KK 7-054 

ASSIGN 715 TO JJ 7-055 

WRITE OUTPUT TAPE 6,716 7-056 

WRITE OUTPUT TAPE 6 , 703 , AVI , AWl , AR 1 , A A1 7-057 

WRITE OUTPUT TAPE 6 , 704 , AV2 , AW2 , AR2 , A A2 7-058 

WRITE OUTPUT TAPE 6, 705, ERR 7-059 

716 F0RMAT(24H FINAL PARAMETER VALUES) 7-060 

WRITE OUTPUT TAPE 6,710 7-061 

IF( KDCAL) 75 1,750, 751 7-062 

750 N=NTT 7-063 

NT1=NTT 7-064 

KGDE7=1 7-065 

GO TO 706 7-066 

751 HH(1)=AT1 7-067 

DO 714 1=2, NTl 7-068 

714 HH( I )=HH( I-l )+DTl 7-069 

N=NT1 7-070 

K0DE7=1 7-071 

GO TO 706 7-072 

715 IF (KDPLT) 733,734,733 7-073 

733 CALL PL0T7(HS, XX, NTT, HH,XN, NTl, KDCAL) 7-074 

734 RETURN 7-075 

END 7-076 



SUBROUTINE PLOT? ( HS »XX *NTT *HH »XN »NT1 >KDCAL) P-000 

C CALLING PROGRAM FOR PLOT ROUTINE P-OOl 

DIMENSION HS( 90) ♦XSC 90) * HH ( 90 ) » XX ( 90 ) , XN ( 90 ) » X ( 18 0 ) » Y ( 1 8 0 ) P-002 

DIMENSION P(11)>K(14) P-003 

C=2. 3025851 P-004 

DO 1 1=1, NTT P-005 

1 XS( I )=-LOGF(XX( I ) )/C P-006 

DO 2 I=1*NT1 P-007 

2 XN( I )=-LOGF(XN( I ) )/C P-008 

C SET UP CONDITIONS TO ENTER PLOT SUBROUTINE P-009 

P{1)=5. P-010 

READ INPUT TAPE 7 , 1 00 , { P ( I ) • I =6 , 1 1 ) P-011 

K(l)=48 P-012 

K(2)=2 P-013 

K(3)=NTT P-014 

K(5)=NTl’ P-015 

DO 3 1=1, NTT P-016 

X(I)=XS(I) P-017 

3 Y( I )=HS(I ) P-018 

DO 4 1=1, NTl P-019 

J=NTT+I P-020 

X(J)=XN{I) P-021 

4 Y( J)=HH( I ) P-022 

IF( KDCAL) 5,6,5 P-023 

5 WRITE OUTPUT TAPE 6,103 P-024 

103 F0RMAT(49HPTPL0T OF ELASTIC CROSS SECTION VS CM ANGLE (DEG), P-025 

1 5X,58HN0TE - EXPERIMENTAL AND CALCULATED ANGLES ARE NOT THE SAME) P-026 

GO TO 7 P-027 

6 WRITE OUTPUT TAPE 6,105 P-028 

105 F0RMAT(49HPTPLDT OF ELASTIC CROSS SECTION VS CM ANGLE (DEG), P-029 

1 5X,54HNDTE - EXPERIMENTAL AND CALCULATED ANGLES ARE THE SAME) P-030 

7 CALL PLDTMY(X,Y,K,P) P-031 

WRITE OUTPUT TAPE 6,104 P-032 

104 F0RMAT(16HPL* EXPER I MENTAL 5X , 14H + CALCULATED) P-033 

100 FORMAT(6F10,0) P-034 

RETURN P-035 

END P-036 
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SUBROUTINE EIGHT 8-000 

C PROGRAM EIGHT SEARCHING PROCEDURE 8-001 

CO'^^'ON AAl ,AA2 ,ACC,ADD,Ak:k:,AK5,ARC»ARM,ARR ,AR1 ,AR2»ASS»AT1 ,AV1 * C-01 

1 AV2 ,AW1 ,AW2,AZD»AZZ,A2D»A3D*A6D » C-02 

2 BB1,BB2»BB3»BB4,BB5»BB6»BB7,BB8*BB9»B10,B11*B12.BM,BP* C-03 

3 CC,CN»CS» C-04 

4 DAI ,DA2 ,DD»DDS,DR1 ,DR2 *DT1 »DV1 *DV2 *DW1 ,DW2 ,D12 * C-05 

5 ERR,FF,GG»HH,HS, C-06 

6 IIl,II2,II3,II4,ICC»INKDDE,IRM*IRRfI3I » C-07 

7 KDCALfKDGRD,KDDLT,KDDE2,<0DE6,KnDE7,'<r!DF8,K800* C-08 

8 LL,LLM,LRM, C-09 

9 NA1X,NA2X»NN1 *NN2,NN3,MN4*NN5,NN6,MN7,NM8 >NNN,NR1X*NR2X C-IO 

COMMON NTN,NTT,NTl,NVlX»MV2X,NV/lX,NW2X*N3NfPL , C-11 

1 QQ,RR ,SN,SS,TT , VC , VV > WW , WX » XC , XN » XX C-12 

DIMENSION BM(51) ,BP(51) , D-Ol 

1 CC( 2,100) ,CN(90) ,CS(51) * D-02 

2 DD( 2,100) ,FF(2,51 ) ,GG(2 ,51 ) ,HH(90) ,HS(90) ,PL( 51) , D“03 

3 QQ(50I ,RR(500) ,SS{500) ,SN(5X) ,TT(500) ,VC(500) ,VV(500) , D-04 

4 WW( 500) ,WX ( 90 ) , D-05 

5 XC(90) ,XN{90) ,XX(90) D-06 

IF(K800) 899,898,899 8-002 

898 K800=l 8-003 

ASSIGN 800 TO III 8-004 

899 GO TO III, (800,811,812,813,814,815,816,817,818) 8-005 

C SET ENTRY CHANNEL AND EXIT CHANNEL 8-006 

800 Z=ERR 8-007 

ASSIGN 890 TO II3 8-008 

IF(DVl) 801,842,801 8-009 

801 ASSIGN 811 TO III 8-010 

ASSIGN 821 TO II2 8-011 

J =0 8-012 

GO TO 811 8-013 

842 IF(DWl) 802,843,802 8-014 

802 ASSIGN 812 TO III 8-015 

ASSIGN 822 TO II2 8-016 

J=0 8-017 

GO TO 812 8-018 

843 IF(DRl) 803,844,803 8-019 

803 ASSIGN 813 TO III 8-020 

ASSIGN 823 TO II2 8-021 

J =0 8-022 

GO TO 813 8-023 

844 IF(DAl) 804,845,804 8-024 

804 ASSIGN 814 TO III 8-025 

ASSIGN 824 TO I 12 8-026 

J=0 8-027 

GO TO 814 8-028 

845 IF(DV2) 805,846,805 8-029 

805 ASSIGN 815 TO III 8-030 

ASSIGN 825 TO II2 8-031 

J=0 8-032 

GO TO 815 8-033 

846 IF(DW2) 806,847,806 8-034 

806 ASSIGN 816 TO III 8-035 

ASSIGN 826 TO II2 8-036 

J=0 8-037 

GO TO 816 8-038 

847 IF(DR2) 807,848,807 8-039 

807 ASSIGN 817 TO III 8-040 

ASSIGN 827 TO II2 8-041 

J=0 8-042 


41 



GO TO 817 8-043 

848 IF(DA2) 808,S49*808 8-044 

808 ASSIGN 818 TO III 8-045 

ASSIGN 828 TO II2 8-046 

J=0 8-047 

GO TO 818 8-048 

C ENTRY CHANNELS 8-049 

811 E=AV1 8-050 

F=DV1 8-051 

GO TO II3»(890,891,892) 8-052 

812 E=AW1 8-053 

F=DW1 8-054 

GO TO II3»(89C,891,892) 8-055 

813 E=AR1 8-056 

F=DR1 8-057 

GO TO II3»(a90,891,892) 6-058 

814 E=AA1 8-059 

F=DA1 8-060 

GO TO TI3»(890, 891,892 ) 8-061 

815 E=AV2 8-062 

F=DV2 8-063 

GO TO I 13, {890,891,892) 8-064 

816 E=AW2 8-065 

F=DW2 8-066 

GO TO 113,(890,891,892) 8-067 

817 E=AR2 8-068 

F=DR2 8-069 

GO TO 113,(890,891,892) 8-070 

818 E=AA2 8-071 

F=DA2 8-072 

GO TO II3,(890,891,892) 8-073 

C EXIT CHANNELS 8-074 

821 AV1=E 8-075 

IF( J) 895,895,842 8-076 

895 continue 8-077 

K0DE8=2 8-078 

RETURN 8-079 

822 AW1=E 8-080 

IF(J) 895,895,843 8-081 

823 AR1=E 8-082 

IF(J) 894,894,844 8-083 

894 continue 8-084 

K0DE8=1 8-085 

RETURN 8-086 

a24 AA1=F 8-087 

IF(J) 894,894,845 8-088 

825 AV2=E 8-089 

IF(J) 895,895,846 8-090 

826 AW2=E 8-091 

IF(J) 895,895,847 8-092 

827 AR2=E 8-093 

IF(J) 894,894,848 8-094 

828 AA2=E 8-095 

IF(J) 894,894,849 8-096 

C FIRST STEP IN SEARCH 8-097 

890 ASSIGN 891 TO 113 8-098 

D = F 8-099 

X = ERR 8-100 

A = E 8-101 

E=E+D 8-102 

GO TO 112,(821,822,823,824,825,826,827,828) 8-103 
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C SECOND STEP IN SEARCH 8-1Q4 

891 ASSIGN 892 TO II3 8-105 

IF(X-ERR) 871»871»R72 8-106 

871 Y=ERR 8-107 

B = E 8-108 

D = -D 8-109 

E = A+D 8-110 

GO TO 1 12 ,(821, 822, 823, 824,825, 826*827, 828) 8-111 

872 Y = X 8-112 

B = A 8-113 

X = ERR 8-114 

A = E 8-115 

E = A+D 8-116 

GO TO 112,(821,822,823,824,825,826,827,828) 8-117 

C THIRD AND SUBSEQUENT STEPS IN SEARCH 8-118 

892 IE (X-ERR) 873,873,874 8-119 

C KEEP GOING IN SAME DIRECTION 8-120 

874 Y=X 8-121 

B=A 8-122 

X=ERR 8-123 

A=E 8-124 

E=A+D 8-125 

GO TO 112,(821,822,823,824,825,826,827,828) 8-126 

C GET BEST ESTIMATE OE E AND ERR AND GO ON TO NEXT VARIABLE 8-127 

873 G=( Y+ERR) /X-2.0 8-128 

IE(G-1.0E-6) 896,896,875 8-129 

875 G=0.5*(Y-ERR)/(X*G) 8-130 

E=A+D*G 8-131 

IE(G) 876,876,877 8-132 

876 ERR=X-(ERR-X)*G*G/(1 .0-G-G) 8-133 

GO TO 896 8-134 

877 ERR=X- ( Y-X)*G*G/ ( 1.0+G+G ) 8-135 

896 ASSIGN 890 TO II3 8-136 

J=1 8-137 

WRITE OUTPUT TAPE 6, 888 , ERR, E 8-138 

88 8 FORMAT ( 1H+27X ,24H INTERPOL AT I ONS , ERRDR = 1 PEI 5 . 8 , 1 3H , PARAMETER= 8-139 
11PE15.8) 8-140 

GO TO II2, (821,822,823,824,825,826,827,828) 8-141 

C END OF A CIRCUIT 8-142 

849 IE( ABSF(X-Z) /X-BIO) 881,880,880 8-143 

880 NN7=NN7-1 8-144 

IE(NN7) 881,881,879 8-145 

879 114=114-1 8-146 

IF (II4) 882,882,800 8-147 

882 II4=NN6 8-148 

DV1 = DV1^^BB9 8-149 

DW1=DW1*8B9 8-150 

DR1=DR1*BB9 8-151 

DA1=DA1*BB9 8-152 

DV2=DV2*BB9 8-153 

DR2=DR2*BB9 8-154 

DA2 = DA2«-BB9 8-155 

DW2=DW2*BB9 8-156 

GO TO 800 8-157 

881 NN5=0 8-158 

K0DE8=1 8-159 

RETURN 8-160 

END 8-161 



APPENDIX B 


SAMPLE DATA SET 
Sample Input Data Listing 


TEST 

CASE# 

CU-ALPHA, 

39.8 MEV* 

SEARCH ON 

R2f A2 

20 

4 

1 10 1 

1 1 

10 40 

35 

1 

1 

1 1 1 

1 1 

1 


48. 

0. 

6.82 

.5 

0. 

0. 

0. 

14. 

6. 

.7 

0. 

0. 

39.8 

6.8 

.1 

9. 

4. 

2. 

20. 

2. 





5. 

OtO 

2.+.0 

l.+O 

O.+O 

1.-7 

. 

5 + 0 

1.-2 

l. + O 



24. 

400. 

1. 

27. 

180. 

1. 

29.5 

77. 

1. 

30.5 

37. 

1. 

32. 

53. 

1. 

35.5 

63. 

1. 

37.5 

40. 

1. 

39.5 

15. 

1. 

41.5 

4. 

1. 

43. 

10. 

1. 

46. 

12. 

1. 

48. 

9. 

1. 

50. 

3. 1 

1. 

52. 

1. 

1. 

54. 

2.5 

1. 

56. 

3.6 

1. 

58. 

3. 1 

1. 

60. 

1.8 

1. 

61.5 

462 

1. 

63.5 

.53 

1. 

66. 

1. 

1. 

68. 

1.1 

1. 

70. 

.8 

1. 

72. 

.45 

1. 

74. 

'.17 

1. 

76. 

.18 

1. 

78.6 

.28 

1. 

79.5 

.22 

1. 

81.6 

.17 

1. 

84. 

.091 

1. 

86. 

.08 

1. 

88. 

.09 

1. 

90. 

• 1 

1. 

91.5 

.095 

1. 

94. 

b05 

1. 




5. 

-30. 

1. 

6. 

0. 

1. 


1 


1 


0 . 

• 1 

63. 


0 . 

-.02 

29. 


1.-8 


1.-12 


5.-4 
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Sample Output Data Listing 


ELASTIC SCATTERING 


TBST 

case;# cu-alpha. 

39.8 HEV 

• SEARCH ON 

R2 

INPUT 

TO PART ONE 





AVI 

X 

4.?999»995E 01 

AW1 

X 

0. 


DVI 

s 

0. 

DW1 

» 

0. 


AV2 

- 

0. 

AW2 

X 

1.39999999E 

01 

OV2 

= 

0« 

DW2 

» 

0. 


AEE 

» 

3.97999993E 01 

ARC 

- 

6.79999995E 

00 

AMP 

s 

4*OOOOOOOOE 00 

A2P 

X 

2.00000000E 

00 

ATI 

s 

1.999999996 01 

OT1 

- 

2.00000000E 

00 

BB1 

s 

5.00000000E 00 

BB2 

s 

2.00000000E 

00 

BBS 

s 

9.999999886-08 

BB6 

X 

9.99999988E-09 

BB9 

s 

5.000000006-01 

BIO 

X 

9.99999988E- 

03 

NN1 

X 

20 

NN2 

X 


4 

NNS 

s 

1 

NN6 

X 


1 

NKODE 

X 

1 

KOPLT 

X 


1 


A2 


AR1 

X 

6.81999999E 00 

AA1 

X 

5.00000000E-01 

DR1 

= 

0. 

0A1 

= 

0. 

AR2 

X 

6.00000000E 00 

AA2 

= 

6.99999994E-01 

DR2 

X 

0.99999999E-01 

0A2 

=-1.99999997E-02 

DEL 

- 

0.99999999E-01 

ARR 

= 

9.000000006 00 

AMT 

- 

6.29999995E 01 

AZT 

X 

2.89999998E 01 

NTl 


40 

NTT 

= 

35 

BB3 

X 

1.00000000E 00 

BB4 

X 

0. 

BB7 

X 

9.99999976E-13 

BBS 

X 

4.99999994E-04 

811 

X 

l.OOOOOOOOE 00 




NN3 

= 

1 

NN4 

X 

10 

NN7 

- 

1 

NNS 

X 

10 

KOGRO 

X 

0 

KDCAL 

X 

1 
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CM ANGLE (DEG) 

24.000 

27.000 

29.500 

30.500 

32.000 

35.500 

37.500 

39.500 

41.500 

43.000 

46.000 

48.000 

50.000 

52.000 

54.000 

56.000 

58.000 

60.000 

61.500 

63.500 
66.000 
68.000 

70.000 

72.000 

74.000 

76.000 

78.500 

79.500 

81.500 

84.000 

86.000 
88.000 

90.000 

91.500 

94.000 


CROSS SECTION (MB) 
4.00000000E 02 
1.80000000E 02 
7.69999993E 01 
3.69999996E 01 
5.29999995E 01 
6.29999995E 01 
3.99999997E 01 
1.49999999E 01 
4.00000000E 00 
0.999999991: 01 
1. 19999999E 01 
9.00000000E 00 
3.09999999E 00 
l.OOOOOOOOE 00 
2.50000000E 00 
3.59999999E 00 
3.09999999E 00 
1.80000000E 00 
6.19999993E-01 
5.29999989E-01 
l.OOOOOOOOE 00 
1.09999999E 00 
7.99999994E-01 
4.49999994E-01 
1.69999999E-01 
1.80000000E-01 
2.79999995E-01 
2. 19999999E-01 
1.69999999E-01 
9.09999990E-02 
7.99999988E-02 
8.99999988E-02 
0.99999999E-01 
9.49999988E-02 
4.99999994E-02 


ERROR WEIGHT FACTOR 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 

l.OOOOOOOOE 

00 
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OUTPUT PROM PROGRAM 


AKK = 2.59434906E 

00 

ADO = 

3.89999992E-02 

AZZ = 2.89627042E 00 

ACC = 4.76877314E-04 

ARR = 9.00899982E 

00 

ARM - 

1. 15829995E 01 

IRR = 231 

IRM = 264 

LRM = 

30 





RECURRENCE RELATIONS 

FOR 

COULOMB 

WAVE FUNCTION 

HAVE FAILED, INCREASE RMAX 



FINAL 

VALUE 

ARM = 1.25189994E 

01 

IRM = 276 

FINAL 

VALUE 

ARC = 6.78599989E 

00 

IRC = 174 


MMM 

St 


0 

LRM = 

30 

LMAX = 28.000E 

00 



AVI 

= 

4-79999995E 

01 

AWl * 0. 


ARl = 6.81999999E 

00 

AA1 

= 5.00000000E-01 

AV2 

= 

0. 


AW2 = 1.39999999E 

01 

AR2 = 6.00000000E 

00 

AA2 

= 6.99999994E-01 

ERROR 

s 

5.804632 13E 

00 







MMM 



0 

LRM = 

30 

LMAX = 28.000E 

00 



AVI 

= 

4.79999995E 

01 

AWl = 0. 


ARl = 6.81999999E 

00 

AAl 

= 5.00000000E-01 

AV2 


0. 


AW2 = 1.39999999E 

01 

AR2 = 6.09999996E 

00 

AA2 

= 6.99999994E-0I 

ERROR 

= 

6.73265034E 

00 







MMM 



0 

LRM = 

30 

LMAX = 28.000E 

00 



AVI 

= 

4.79999995E 

01 

AWl = 0. 


ARl = 6.81999999E 

00 

AAl 

= 5.00000000E-01 

AV2 

s 

0. 


AW2 = 1.39999999E 

01 

AR2 = 5.89999998E 

00 

AA2 

= 6.99999994E-01 

ERROR 

= 

5.47076815E 

00 







MMM 



0 

LRM = 

30 

LMAX = 28.000E 

00 



AVI 

= 

4.79999995E 

01 

AWl = 0. 


ARl = 6.81999999E 

00 

AAl 

= 5.00000000E-01 

AV2 


0. 


AW2 = 1.39999999E 

01 

AR2 = 5.79999995E 

00 

AA2 

= 6.99999994E-01 

ERROR 

= 

5.55709577E 

00 

INTERPOLATIONS, 

ERROR= 5.45254004E 00, 

PARAMETER= 

= 5.87054473E 00 

MMM 

sc 


0 

LRM = 

30 

LMAX = 28.000E 

00 



AVI 

= 

4.79999995E 

01 

AWl = 0. 


ARl = 6.81999999E 

00 

AAl 

= 5.00000000E-01 

AV2 


0. 


AW2 = 1.39999999E 

01 

AR2 = 5.87054473E 

00 

AA2 

= 6.79999989E-01 

ERROR 

* 

5.08810759E 

00 







MMM 

X 


0 

LRM = 

30 

LMAX = 28.000E 

00 



AVI 

X 

4.79999995E 

01 

AWl = 0. 


ARl = 6.81999999E 

00 

AAl 

= 5.00000000E-01 

AV2 

= 

0. 


AW2 = 1.39999999E 

01 

AR2 = 5.87054473E 

00 

AA2 

= 6.59999985E-01 

ERROR 

= 

5.02560300E 

00 







MMM 

— 


0 

LRM = 

30 

LMAX = 28.000E 

00 



AVI 

= 

4.79999995E 

01 

AWl = 0. 


ARl = 6.81999999E 

00 

AAl 

= 5.00000000E-01 

AV2 

X 

0. 


AW2 = 1.39999999E 

01 

AR2 = 5.87054473E 

00 

AA2 

= 6.39999986E-01 

ERROR 

= 

5.28400493E 

00 

INTERPOLATIONS, 

ERROR= 5.01065475E 00, 

PARAMETER^ 

= 6.66104484E-0I 

MMM 



0 

LRM * 

30 

LMAX = 28.000E 

00 



FINAL 

PARAMETER VALUES 







AVI 


4.79999995E 

01 

AWl = 0. 


ARl = 6.81999999E 

00 

AAl 

= 5.00000000E-01 

AV2 


0. 


AW2 = 1.39999999E 

01 

AR2 = 5.87054473E 

00 

AA2 

= 6.66104484E-01 


ERROR = 5.01065475E 00 


4 ? 



/ 


CALCULATED VALUES OF ELASTIC CROSS SECTIONS 


CN ANGLE (OEG) 
20*000 
22.000 

24.000 

26.000 

28.000 

30.000 

32.000 

34.000 

36.000 

38.000 

40.000 

42.000 

44.000 

46.000 

48.000 

50.000 

52.000 

54.000 
56*.000 

58.000 

60.000 
624000 

64.000 

66.000 

68.000 

70.000 

72.000 

74.000 

76.000 

78.000 

80.000 
82 w 000 

84.000 

86.000 

88.000 

901.000 

92.000 

94.000 

96.000 

98.000 


EL XSCTN (MB) 
9.62958920E 02 
5.70692205E 02 
4.88)53476b 02 
3.47278482E 02 
1.54765056E 02 
5.00323659E 01 
5.44254661E 01 
B. I8312800E 01 
6.98226541E 01 
3. 14305583E 01 
6.26571721E 00 
6.38151097E 00 
1.50700368E 0) 
1.59072275E 01 
8.70288706E 00 
2.1006)249E 00 
8. 18222284E-01 
2.62395263E 00 
3.55630508E 00 
2.49202803E 00 
9.052n937E-01 
2.32998583E-01 
4.70808333E-01 
7.95787382E-01 
7.09919691E-01 
3.67436954E-01 
1.28562714E-01 
1. 12168066E-01 
1.82845095E-01 
I.93067770E-0T 
1.29695812E-01 
6.34438616E-02 
3.98642811E-02 
4.57365376E-02 
4.89706349E-02 
3.83481470E-02 
2.40016937E-02 
1.66017307E-02 
1.58055563E-02 
1.52522676E-02 


COUL XSCTN (MB) 


3.42673385E 

03 

2.35053077E 

03 

1.66742258E 

03 

1.21677208E 

03 

9.096215I3E 

02 

6.94346786E 

02 

5.39770323E 

02 

4.26402742E 

02 

3.41689688E 

02 

2.77327996E 

02 

2.27695602E 

02 

1.88905549E 

02 

1.582)94026 

02 

1.33674738E 

02 

1. I3843484E 

02 

9.76716113E 

01 

8.43710697E 

01 

7.33456379E 

01 

6.4I393006E 

01 

5.639965426 

01 

4.98518568E 

01 

4.42796361E 

01 

3. 95114)036 

01 

3.541003026 

01 

3.18650621E 

01 

2.87869832E 

01 

2.610277096 

01 

2.375249776 

01 

2.16866997E 

01 

1.986434 I2E 

01 

1.8251 1987E 

01 

1.68185945E 

01 

1.55423896E 

01 

1.4402 1747E 

01 

1.338062156 

01 

1.24629647E 

01 

1.16365713E 

01 

1.08905929E 

01 

1.021568606 

01 

9.603774196 

00 


EL/COUL 
2.81013629E-01 
2.42792907E-0) 
2.92759302E-0I 
2.85409638E-0I 
1.70142256E-01 
7.20567405E-02 
1.00830786E-01 
1.91910774E-01 
2.04345217E-01 
1. 13333519E-01 
2.75179541E-02 
3.37814897E-02 
9.52477157E-02 
1. 18999504E-01 
7.64460701E-02 
2. 15068883E-02 
9.69790089E-03 
3.57751691E-02 
5.54465818E-02 
4.41851646E-02 
1.81580383E-02 
5.26198047E-03 
1. 19157560E-02 
2.24735016E-02 
2.22789359E-02 
1.27639964E-02 
4.92525154E-03 
4.72236919E-03 
8.43120873E-03 
9.7193I386E-03 
7.106153016-03 
3.77224508E-03 
2.56487462E-03 
3.17566884E-03 
3.65981761E-03 
3.07696822E-03 
2.06260872E-03 
1.52441016E-03 
1.54718500E-03 
1.58815350E-03 
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Sample Output Data Plot 


PLOT OP ELASTIC CROSS SECTION VS CM ANGLE (DEG) NOTE - EXPERIMENTAL AND CALCULATED ANGLES ARE NOT THE SAME 

0 10 . 20 . 30 . 40 . 50 . 60 * 70 . 60 . 90 . 100 . 

3.0 


2.0 


1.0 


* » 

1 * • 

+ 


+ ♦ * 


1.0 


♦ 1 


» — 


2.0 


10 . 20 . 


30 . 


40 * 


50 . 


60 . 


70 . 


80 . 


90 . 


100 . 


CO 


• EXPERIMENTAL 


♦ CALCULATED 
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NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 
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